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ABSTRACT 


Structure in the conductance o of normal metal- 
insulator-metal junctions at very low bias is explained 
through a nonequilibrium treatment of the tunneling 
process. In particular, the related peak in the deri- 
vative do/dV is quantitatively accounted for by the 
blocking of otherwise-available electron tunneling 
states due to the finite electron relaxation rates in 
the metal electrodes. A transport model which includes 
the effects of elastic and inelastic scattering of 
electrons in the metal electrodes is used to obtain 
expressions for o and do/dV. The results are compared 
with the temperature dependence of the peaks of do/dv 
obtained from measurements on Al1-Al and Pb-Au tunnel 
junctions. The nonequilibrium theory is extended to 
the case of tunneling in superconducting junctions and 
results in a non-discontinuous rise in the tunneling 


current at the gap edge. 
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CHAPTER I 


INTRODUCTION 


1.1 General Background 


A transport phenomenon which has been of particu- 
lar interest over the past decade is the transport of 
electrons through thin insulating films. One system 
(commonly called a tunnel junction) that is used to 
study such effects is made by depositing a metal film 
onto a glass siirde, Jetting it oxidize for a few 
Minutes and then depositing a second film crossing the 
first. Such an oxide film on metals like aluminum, lead 
or tin will ordinarily be continuous and may be of the 
order of 20 A thick. The two metal films are therefore 
not in electrical contact; yet if a voltage is applied 
across the two, a current is found to flow and is in 
fact proportional to the applied voltage. This is the 
behavior expected if the current was conducted through 
small metal bridges in the oxide. However, when the 
metallic films are superconducting, it becomes clear 
that the current is carried by a tunneling mechanism. 
To understand the phenomenon of electron tunneling through 
thin insulating films, we must understand the behavior 
of an electron in the vicinity of a metal surface. We 


shall return to tunneling in superconductors in Chap. V. 
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Theories of metal surfaces adjacent to a vacuum 
or an insulator have, relatively speaking, lagged far 
behind the bulk theories of metals which are now capa- 
ble sormgivying GUuantitatively accurate descriptions of 
Wldemchasses Of metalsy This) is primarily due to the 
loss of translational invariance at the surface of the 
metal and the rapid decrease of electron density there. 
In contrast, the translational invariance of the lattice 
in the bulk metal introduces important elements of sim- 
DIletey into the calculatvons: of the bulk cheorist: 

The earliest model calculations of the metal inter- 
face were made by Frenkel (1928). He replaced the 
lattice of the metal by a constant positive "jellium" 
density which ended as a step discontinuity at the metal 
surface. He applied the Thomas-Fermi Method together with 
Poisson's Equation to calculate the electron density 
near the interface. Integrated charge neutrality as well 
as the vanishing of electronic charge density at large 
distances exterior to the metal were used as boundary 
Condtiteons in has calculations. Bardeen (1936) used the 
same approach except that he used the many-electron 
Schrodinger equation rather than the Thomas-Fermi Appro- 
ximation to relate the electron potential to its charge 
density. Recently, a new formulation was developed in a 
series of papers by Kohn and his collaborators (1964-1971). 


Kohn uses the (presumed known) results for the ground 
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state energy of a uniform system to calculate the energy 
of an inhomogeneous system by writing the latter as a 
Punctironal vor the Local electron density. 

The picture of the metallic interface that has 
emerged from these studies is very similar to what one 
feels should happen intuitively. The potential an 
electron sees inside a metal is periodic with the 
period of the lattice, going through deep depressions 
in the neighborhood of each atomic nucleus. It rises 
sharply at the surface, forming a potential barrier 
which keeps the electrons confined to the metal. This 
transition zone at the surface is of the order of an 
electron Fermi wavelength wide. The potential energy 
within a metal is lowered due to the actions of the 
attractive restoring forces of the atomic nuclei and 
the exchange forces between electrons of the same spin. 
thes translational invariance of the bulk, lattice forces 
the electrons in a solid to be grouped into energy bands 
separated by gaps of forbidden energy. The lower energy 
bands are filled with electrons, each energy level in a 
band having its full quota of two electrons according to 
the Pauli Exclusion Principle. In metals the highest 
energy band is only partially filled and forms the con- 
duction band, only electrons of this band participating 


in conduction and transport processes. 
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If we consider a conduction band electron wave- 
Function in one of the metal strips of a tunnel Junction, 
it is clear that the wavefunction does not drop immed- 
iately to zero at the interface with the oxide. Rather 
it will penetrate and decay exponentially into the oxide. 
The wavefunction may still have a very small but finite 
value at the opposite edge of the oxide and therefore it 
1s possible, for the. electron. to ,makesa, transition, or to 
"tunnel" into the second strip. This is the same behavior 
that would be anticipated if the oxide were to be replaced 
by an equivalent thickness of vacuum. The conduction 
electrons of interest would have energy below the vacuum 
level, and therefore the wavefunction would have negative 
kinetic energy within the vacuum and again would decay 
exponentially. This is precisely the situation described 
as quantum-mechanical tunneling and treated in elementary 
quantum=mechanics textbooks. (e.g. see Schifft for, instance) . 
Oxides are predominantly used as the barriers in tunnel 
junctions because of their ease of fabrication. In order 
to easily observe a tunneling current between two metals 
separated by a vacuum, the metals should, in practice, 
be spaced less than 50 A apart .«, The, construction efpsuch 
a tunnel junction presents a very difficult task since 
the smoothness of ordinary metal surfaces is certainly 


not of this order, even when they are specially prepared. 
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The only difference with the oxide layer is that the 
exponential decay of the electron arises from the fact 
that the energies of interest lie within the forbidden 
energy band of the oxide rather than from negative 


kinetic energy as in vacuum. 


1.2 Historical Resumé of Tunneling Investigations 


The first application of the idea of quantum- 
mechanical tunneling to describe emission from a bulk 
metal was the calculation of Fowler and Nordheim (1928) 
that explained field emission from free electron metals. 
Frenkel (1930) soon extended this calculation to the 
case of current flowing across a voltage biased metal- 
vacuum-metal system. He accounted both for the change 
in character of the barrier potential (which he assumed 
was trapezoidal in shape) and for the thermal occupation 
probabilities of the electronic states in the metal 
electrodes. An extension of Frenkel's calculation was 
made by Sommerfeld and Bethe (1933) to the case in which 
the barrier between the two metal electrodes is an 
insulator rather than a vacuum. The band structure 
and dielectric constant of the barrier region was taken 
to be characteristic of a bulk insulator. 

The quantum-mechanical tunnel current flowing 


through a thin insulating layer between normal metal 
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electrodes was first experimentally investigated in 
detail by Fisher and Giaever (1961). They measured the 
current-voltage characteristics of junctions made with 
evaporated electrodes of aluminum and a thermally grown 
oxide layer. Detarled comparison with the theory of 
Holm (1951) verified that tunneling was in fact taking 
place. When their devices were subsequently cooled to 
superconducting temperatures it was found that super- 
conducting energy gaps and densities of quasi-particle 
states characteristic of the bulk films could be obtained 
from the now non-linear current voltage characteristics 
(Giaever, 1960). 

Giaever's work immediately stimulated a sequence 
of measurements on both superconducting and normal metal- 
insulator-metal junctions. The early experiments on 
superconducting tunnel junctions provided additional 
confirmation of the Bardeen-Cooper-Schrieffer Theory of 
SUupeLrconductivity, (1957, herearter referred “co as the 
BCS theory) for weak coupling superconductors (e.g. Al, 
In, Sn). Subsequent experiments also verified the theory 
Of Strong coupling superconductors (e.g. Pb, Hg): The 
history of developments of measurements on superconduct-— 
ing tunnel junctions has been reviewed in detail by 


McMillan and Rowell (1969). 
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Giaever's experiments also provided the motivation 
fon a ,ecxaminattonsotatheatheory,ofstunneling! 1neorder 
to incorporate into it a description of the many body 
properties of the electrodes and thus provide an adequate 
base for theoretical calculations. The models of tunnel- 
ing (and their experimental manifestations) prior to 
Giaever's pioneering measurements had the common feature 
that they were based upon a one-electron description of 
the tunneling process. Although the theory of bulk 
superconductivity had already been developed by BCS, 
the superconducting tunnel conductance, o = dI/dV, seemed 
to be a direct measure of the many body bulk density of 
states in contrast to the predictions of a one electron 
Eheory: 

The subsequent reexamination of tunneling theory 
led to the transfer-Hamiltonian formulation of “the tun- 
neling current by Bardeen (1961) and Cohen, Falicov and 
Phillips (1962). This reformulation of tunneling theory 
led to the prediction by Josephson (1962) and subsequent 
observation by Anderson and Rowell (1963) of phase 
coherence between two superconductors separated by a 
thin tunnel barrier (now commonly known as Josephson 


tunneling). 


_ Oo 7 
-* =. a 
; * 4 ag PT 
: . 7 
aa eat ae oele adiomitegee e'tavesto - 


tshso ni pailersus to’ yrosdt sft to foideniaaxset | 


sat315aed 5 +4 4ek adatoq2oont 
4 afd to sottragona 


_ 
: 4. yhed ynom silt to nol 


- ——“(: ot ePSbs Me SPhivomg Bert ins esbortosis 
> 


7 <fenaus to elabom ent stotjatioisa Lsotsetosts toz cand 


. 
oF robiq (amottst=ast fash Iganambtagqxe tien Bie) ead 


- 

: eivtest nommoo ott bid atcsmsrs59m pret (S504 a! soveeho 
_ 

a 


25 gebiqiistsh aoxdosio-sio & nogy heasd sxew yors Snel 


died tq yroors oie HopodsttA .8250 


,208 yd bsqelsavsb nesd ybsoits bea ya lvatoubaeat 


hatiese .VE\ID = 0 , sonmsyousnos Tontnua po td oubagosogue: of 


7 
36 ydredob alud ybod yasm eit to sivesem Joe2Lh ote 
Deysnes Gl estete- 


moswsosis sie 6 ZO vnoisorhbexq ons .os)4 


— 
ee 


yioods prilenns to foltsiimsxess tretipevedie oat” 
mut gis to noitslunzo2 ns rnotiimesh-rstensis 


bas yoourlse? (aede> bas (eel) nesbasa ya Jo9rtve @ ife 
2 
yroe enilennus 46 notdeiumrotss efdt . (Self) aqiil 


7 
-iaaeel edad Bas. (SdCL) fmoerqesot Yo dortpibatg git of 7 


: , 
- o sencig ay (&act) Lfowor bas poerSebart yd noisavTes 
=i 
7 | 5 yd bodaisgse i kai aA ws, asewdod Sone 
16> i 
a omeqonot ee swond acini worn) asizxsd tsanva-i 


1.3 xperimental Procedures in Tunneling 

Since experimental data is presented in this thesis 
and compared with theory, a brief description of the 
measurements and data is given to provide an orientation 
for subsequent theoretical analysis. Unless otherwise 
noted, the experimental data in this thesis is presented 
with the kind permission of Dr. J.G. Adler. 

The stepseiim the? fabrication of a tunnel junetion 
are outlined in esq. (1-1). A, smoothly polished glass 
slide with indzum contacts (a) is) placed in a high 
vacuum environment. A metallic strip about 1 mm wide 
and 1000 to 3000 A thick is deposited across the indium 
contacts (b). The metallic strip is then allowed to 
oxidize to form an insulating layer some tens of Anstroms 
thick (c). A ssecond strip of metal‘is then,evaporated in 
vacuum across the first. one (d). The connections of 
current and voltage leads are shown. The junction is 
then taken out of the evaporator and placed into a liquid 
He* cryostat where measurements down to 1°K can be made. 
A superconducting magnet surrounding the junction may be 
used to subject the junction to magnetic fields as high 
as 60 kilogauss. Normally several films are evaporated 
across the oxidized base film to simultaneously provide 
several junctions for measurement. 

In a typical normal metal-insulator-metal (M-I-M) 


junction measurement, a constant current I, is passed 
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through the junction together with a small modulation 
Currents cossati where <ieispaysmall: constant current 
modulation amplitude, w is 27 times the modulating 
Prequencysandut wsncthettime-2 thes cusrented = It 6cos wt 
will result in a voltage response V(I) across the junc- 
tion. Since the junctions are weakly non-linear, the 


voltage response may be expanded in powers of small 6: 


2 
V(I) = viz.) + ()_ 6 cos wt + $(2Y) 5? cos7wt 
ar’ T 
oO 
“jan Meh sen 6 
av 1,d°v 2 
= V(I,) + (G4) 6 cos wt +<(— 87 (1+ cos wt) 
I fat’ I 
oO 
= ee neBe cies 


Vane Lsetne G.c. bias Of sLhie tunnel junction and 12t tel-— 
lows from the above equation that if 6 is small compared 
to I, and constant, then the component of voltage across 
the junction at angular frequency w is proportional to 
(voir and the component at 2w is proportional to 
(d°v/ar?) | . Both these two components are measured by 
means of eee sensitive instrumentation. The quantities 
Of interest, the current I(V), the dynamic conductance 
o(V) = dI(V)/dv, and derivative of conductance do(V)/dV = 
d*1(v) /av? are then reconstructed numerically from the 


above measurements (see Adler et.al. 1971) and form the 


finaledata ton theoretical interpretation. 
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The objective of the present investigation was to 
carry out an intensive interpretation and theoretical 
analysis of the conductance minima centered about zero 
bias (hereafter called zero-bias anomalies or ZBA) 
occurring in normal M-I-M junctions under experimental 
investigation by Dr. J.G. Adler. 

In Chapter II of this thesis we present a review 
of conventional tunneling theory both from the one- 
electron picture and the perturbation (transfer- 
Hamiltonian) point of view. This review will provide 
the concepts in tunneling required to explain the phy- 
Sical origin of the zero-bias-anomalies. It will also 
provide us with the means for isolating the zero-bias- 
anomalies from the gross background features of electron 
tunneling which tend to obscure these anomalies. This 
will allow us to compare theory with experiment in an 
effective manner. 

The theory of the ZBA as an interference or "block- 
ing" effect in the transport of electrons across the 
junction is presented in Chapter III. We compare theory 
to experiment in detail for two cases; an aluminum- 
aluminum junction and a lead-gold junction. Partial 
results of Chapter III have been published in Physical 


Review Letters (P.N. Trofimenkoff, H.J. Kreuzer, W.J. 
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Wattamaniuk and J.G. Adler, 1972). 

Chapter IV is presented for completeness and 
involves a discussion and phenomenological explanation 
of the nonlinearities that occur in the tunneling process 
due to the presence of the barrier oxides. 

Me-Eicory-ot-quasa particle tunneling in super— 
conducting tunnel junction is reviewed in Chapter V. 

The extension of ideas put forward to explain the ZBA's 
in normal metal junctions leads to a plausible explana- 
tion of the nonideal behavior in the current-voltage 
characteristics of a superconducting-superconducting 
tunnel junction. 

In the final Chapter (VI), we give a general dis- 
cussion of the geometry of the tunnel junctions based on 
related experiments made by Dr. J.G. Adler. It is shown 
that rtunneling occurs at restricted areas, of the tunnel 


Functions: ofeinterest jin «this: thesis: 


> 


i a ze <a 
ferek eprstveantl inant elt 2 7 
biog ne. wis i060 eonessid siz ots 


stoque al omiiviads elo bi isq-Eaetip 30 mens oe? . 
: Vo tadqed? “at bowsives ei Adidonut Lodanst pitt soubH0® 


wv. —" 


2’ Ag® ist sed ade: ad biuswsot sug aie? 23 adianehan z 
~aitel gnS aidievaly 5B ot absel aro itoayt Lasem Lamzos ad 
Sperioy-sasiiw sdy af toiverisd taebinon sit Jo mori 
pitivovbnoSséquerpitsoubaooieqna & 46 axiteissdosmsio : 
| -poisomp os 
-—eib Satis & evip ow (LV) cetqsdd Lenii aa4 of 
no ‘bseed enor7: site Leong ait to yasomosp oils io pie 
| awORe et 41 sa8tbe he Si is Gamt eo ee bSaa | . 
_ ‘techn Si io igbate pedoixsves ts ‘e2up3. oii lenood +2 oan 


laleadt 21nd oi tesposnt Fo me Lteaut 


is 
CHAPTER sl. 


CONVENTIONAL TUNNELING THEORY 


(NORMAL METALS) 


2.1 The One-Electron Picture of Tunneling 


In this section, we derive the tunneling current 
across a M-I-M junction using essentially the arguments 
of Sommerfeld and Bethe (1933) and Fan (1942). It is 
assumed that the insulator may be replaced by an 
effective one-electron potential barrier whose thickness 
is equivalent to that of the oxide. Any distortion of 
band structure in the metals due to mismatch of metal- 
insulator work functions is neglected. Some justification 
for these assumptions has been given in Section 1.1 of 
the Introguction., | ine tunnel junction atselr will be 
envisioned geometrically as an insulator of even thick- 
ness sandwiched between two parallel metals. The metals 
are treated in the free electron approximation in which 
the periodic potential of the crystal lattice is neglected. 
In accordance with this model, only the electrons in the 
ae e filled conduction band may tunnel. A convention 
that is adopted throughout the remainder of this thesis 
is. that quantities characteristic. of the lett and right 
hand metals shall be identified by subscripts & and r 


respectively. 
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The energy~space diagram of two metals M, and M, 


& 
separated by an insulator of thickness s and gap energy 

oe is represented schematically by Fig. (2.1). At 0°K 

all electronic states in the metal are occupied up to a 
-critical value of the kinetic energy u (the Fermi energy). 
If metals My and M are different and are shorted out by 
an external conducting wire they will reach (by electron 
exchange) an equilibrium state in which the chemical po- 
tentials in both metals are equal (i.e. the Fermi levels 
coincide). In this state an intrinsic voltage ve (the 
contact potential) arises so that eVvo= ee where e 

is the electron charge and x is the metal-insulator work 
junction (i.e. the work required to displace an electron 
at the Fermi level in the metal into the conduction band 
of the insulator). We shall measure all energies from the 
bottom of the conduction band in Me: 

When an external bias V is applied to the junction 
so that M, is the cathode and M the anode, the Fermi 
Level of M,. will be lowered by an amount eV with respect 
EO2tiatleor M,- The change in electron potential energy 
in a direction x normal to the junction is represented by 
the trapezoidal barrier.(curve-(a)-of- Pigs (2sl))\ 2i-the 
interaction of the tunneling electron with the metals is 
neglected. However, the abrupt change in the potential 


barrier is only an ideal model and is physically unreali- 


zable since abrupt changes in potential imply infinite 
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Fields. In fact, we expect that the barrier potential 
changes smoothly (curve(b), Fig.(2.1)) because the 
transition region between the metal and insulator must 
be of finite extent. The electron will also experience 
an attractive image force in the barrier region due to 
the metal surface (Simmons, 1969). The voltages of 
interest in this thesis are of the order of millivolts 
as contrasted to the Fermi energies and potential barrier 
heights which are of the order of volts. In addition 
the low temperatures of interest here ensure us that 
the Fermi surfaces are sharp. Thus electron transfer 
across the barrier will occur entirely by tunneling at 
Fermi energies. 

The electron current density incident on the metal- 


insulator interface from the left is given by 
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where V is the volume and fk is the momentum of an elec- 
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m is the free electron mass and f° is the equilibrium 


Fermi-Dirac distribution function defined by. 


oe) » ak 
(Ey = yeas ’ (2.4) 


l1l+e 


where kp is Boltzmann's constant and T is the temperature 
im “K. The tactor, of Zein front. of the stmmation over 
free electron states accounts for spin while ve is the 
positive x-directed component of electron velocity 

ie = Ak, /m) . With a large number of states and box 


quantization 


k on 
therefore 
w2e Ue O 
Ue ae I dk | | dk l=) i oa a fa ev) P 


(2:5) 


where dk = Gioia. 


Now let the wave vector k of an electron in the 


left iano metalschangea, to kK “as it -cConnels into the 


right hand metal. We shall assume that energy is con- 
served (1.e. Ey = Evo) along with the parallel component 
of momentum (i.e. k, = k' Le 


8 aco 
The latter simply corresponds to specular reftlec- 


tion at the metal-insulator interface. Now whether or 
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nocean electron tunnels into the right hand metal 
lepenas=upolr: 
(i) the availability of an empty state of energy 


E, in the right hand metal. This is speci- 


k 


fied by the distribution of holes there, 


oO 
a (Ey - HL) 7 


(ii) the probability of the electron penetrating 
through the barrier into the right hand side. 
We specify the probability by a transmission 
coetrrerent Pk) whieh will obveotrse-depend 
upon the bias and explicit shape of the barrier 


potential. 


The transmitted current density from left to right is 


thus given by 


2e use re) 
af = 3| ok | | dk (aie i (Ey) -u ev) 
(217) k_ >0 a 
So 


9 eae) 


O 
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where E, is given by Eqn. (2.2). 
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Similarly, there will be a transmitted current from 
the right hand metal into the left hand metal which by 


(20) > | k'>0 
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where now Ey = A7k'7/2m (remember that we are measur- 
Ing enerqiesetrom the bottom of «the sazght hand band) . 
The above equation for the backcurrent must be modified 
slightly since a’ glance at Fig. (2.1) -shows that “an 
electron in the left hand metal cannot tunnel into a 
right hand state that is below the bottom of the band 
on the right hand side. In fact, it is easy to show 
that in such a case, the energy and momentum conservation 
relations that we have imposed lead to an imaginary x 
component of the wave vector there, implying that the 
wavefunction decays there. This is of course to be 


expected and may be rectified by imposing the following 


condition 

Pes; ky) =i) 
ie 

ae See ue ati. + eV) . 
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where 


Bone. (2.0), NOW Ver similar structurally tovmon.. (26) 5 


At zero bias (eV = 0) no current can flow, therefore 


we must have 


This implies that 


)=Pt (fae + BG -utev) » ok). (2.9) 


PAK i 


If we now combine the forward current (Eqn. (2.6)) with 
the backward current (Eq. (2.8)), the net current from 


left to right is given by 


Hk 
we 2. fa | dk, (—) [f (E, =u, ev) - cP Hiei) ]P(k), 
k >0 = a 
x (20) 
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translationally invariant in the parallel plane (i.e. 


If the barrier is 


where Ey =u 


nS Hyt ev + 
ya Diane) Chen the transmission, coer ficient, Pk) ucan 

only depend on ky or the x directed kinetic energy 

E. = fi*k2/2m available for tunneling. Let us in addition 
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define the parallel kinetic energy of the electron by 
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Thus 
Ak, 27m i { 
{ax, | Se armel, re, | dE | aa ; 
>0 0 0 


and the tunneling current density now becomes 
fe) fe) 

| P(E) dE, | GEy Ee (BL TE Hy) £ (E+E 

0 0 


47em 
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(20h)? 


pH yg tev) I. 


Caters) 


The integral over E,, may easily be evaluated using Egn. 


|| 
(2.4). The result is commonly called the "supply func- 


tron. N(EL), ise: 


(e) 
ie cle (E+E yu, tev) J 


= (©) 
N(E,) = | ey (E+E 
0 


Lerheexe (ieee) 
= k,T in | ——_—4_= 3 _ = pel Bee 
lot exp (u,-E,-eV/k,T) 


Thus the current density now assumes the simple form 


_ 47em 
J = | P(EL)N(E,) dE, j (2043) 
0 


where h = 27h is Planck's constant. The supply function 
N(E,) simply represents the difference between the number 
of electrons per unit time and area, incident on opposite 
sides of the barrier, and having x-directed kinetic energy 


in the range E, to E,+dE,- It is entirely statistical in 
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nature and restricts the tunneling process to those elec- 
trons whose total kinetic energy is close to the Fermi 
energy. The breakdown of the tunneling current into an 
integral over the product of geometrical and statistical 
factors is very convenient and provides the starting 
point for any quantitative analyses of the current- 
voltage characteristics in tunnel junctions. If we 
define the total kinetic energy by E = EL+ Evy then 


Edt. (2eLl). May. oe: Casteincomthie, Lorn 


=} 


j= Stem | dE [£"(E-y,) - £~ (E-y,+eV) ] | P(E) GE es. (2. 1A) 
h 
0 0 


This equation displays more clearly the fact that tunnel- 
ing is restricted to energies close to the Fermi energy. 


At zero temperature we have 


| 
be 


£° (E-uy) = £°(E-u, +eV) PSS) ee ey 


0 otherwise . (2045.) 


Therefore Eqn. (2.14) becomes 


a 
T(V) |peg = =a (eV) | P(E,) dE, (2.16) 
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provided |lev| << Hp + 


For small bias we expect that P(E, ) will be fairly in- 
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sensitive to the voltage. The tunneling current expressed 
by sHon. (2.16) thengshows that the -tunnel junctzon is 


Ohmic for very small bias as confirmed by experiment. 


2.2 Barrier Models and Transmission Coefficients 


We now proceed to compute and discuss the probabi- 
lity of an electron tunneling through the oxide layer. 
The most natural approach is to construct incident and 
reflected waves in the left hand metal, match them to 
exponentially decaying waves of the same energy in the 
oxide, and finally to match these to a transmitted wave 
in the right hand metal. 

We shall consider two different models of the 
potential barries,_the first of which is illustrated in 
Fig. (2.2a) and is characterized by sharp walls with a 
smoothly varying potential inside. Since we have trans- 
lational invariance parallel to the barrier, we need 
only consider an electron incident on the barrier with 
an x-component of kinetic energy Bx = nok, /2m and 
which emerges on the other side with x-directed kinetic 
energy E,. = figks Jeme The ancident and reflected waves 
will be De the ee a eg while the transmitted 
wave is given by exp (ik, x). The potential inside the 

ca 


barrier is left arbitrary and is denoted by 9$(x). Using 


the WKB approximation (e.g. see Davidov (1965)), the 
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Fig. (2.2) Models of tunneling barrler. 
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wave functions within the barrier are 


x 
exp (+| R(y) dy), 9 where espe ( rea) 
aca) } | ee es 


This approximation is justified provided that 9 (x) 
varies Sufticiently ‘smoothly within the barrier. “The 
conditions of smoothness on the potential may be ob- 
tained by simply substituting the wave function into 
the Schrodinger equation for the electron in the 


barrier. These conditions are simply 


2 
ah dk lL 6d 
less | mat and = —| 
K K ax 


If we proceed to match the wave functions and their 
derivatives at the walls of the barrier, we find that 
the transmission coefficient (the ratio of transmitted 


to ineident particle flux) is given by the formula 
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Consider the case where the barrier potential is constant 
(1,.@. 96(x) = bt S constant) and the bias is small. If 
the metals are similar so that iy See Kon’. mae?) 
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2m 
a E See (> -E,) 


5 a ae) e : (2.18) 


O e) 


P(E,) = Ks) (l 


which is just the probability for a particle to tunnel 
through a square barrier of height %% and width s 
(Davidov, 1965). The transmission probability depends 
exponentially on the width of the barrier as expected. 

A more physical model of the oxide barrier with 
smooth walls ts depicted in’ Fig. (2.2b)... Unfortunately 
it is difficult to deal with such a model rigorously 
and therefore we again resort to the WKB approximation 
under the assumption that the potential is smooth enough. 
The WKB solutions are matched at the turning points Xo 
and X,- The transmission coefficient for such a case 
is given by Kemble (1937) and is discussed in detail by 
Miller and Good (1953). The result is simply 

-2 fr K (x) dx 


P(E.) =e . (2.19) 
x 


where Xo and x, are the turning points (i.e. K (x) ) =K (x) = 


0) and «(x) was defined in the previous model (K(x) = 


[z= (> (x) -E,) ee Because’ of its’ ‘simplicity and ease of 


application, this approximation is frequently used in 


quantitative calculations of the tunneling current. 
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Although both models agree in their exponential 
dependence of the tunneling probability on the thick- 
ness and height of the oxide barrier, they differ 
with respect to their preexponential factors. The 
preexponential factors present in the model with sharp 
walis contain, in principle, information about the band 
structure andsdensity of states of both metals. 
Harrison (1961) has calculated the transmission co- 
efficient for. a Square barrier taking into account the 
Structure of the metal lattice by also matching the 
Bloch parts of the electron wavefunction at the inter- 
face. He concludes that details of the lattice structure 
are also contained in the preexponential factors. How- 
ever, experiments to date have not managed to confirm 
these details. The dependence of the tunneling probability 
on the details of the wavefunctions seems rather unrealis-— 
tic when the boundaries are not ideal as one suspects in 


the case in tunneling junctions. 


Zio Magnitude and Angular Dependence of the Transmission 


Coefficient 


In this section we shall estimate the magnitude of 
the transmission coefficient as well as its dependence 
upon the angle of incidence of the tunneling electron. 


Since the experiments are done at very low temperatures, 
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we shall make use of the current density given by Eqn. 
(2.16) for small bias and 0°K to relate the transmission 
coefficient to the resistance of the junction. We con- 
sider a typical junction with aluminum electrodes and a 
resistance of 50 ohms. From Eqn. (2.16) the tunneling 


current I is given by 


2 


4te MAP .u 
I = ae Ne (2.20) 
h 
where A is the effective tunneling area and Po isa 
dimensionless quantity defined by 
u 
et 
Dare Ti | P(E.) dE, 3 (2 seek) 
0 


LB hes tesistancerote the tunnel junction 2s—_R, <cohen, Lom 


Bomen2 sc), LesLOLLows. that 
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fe ee eee Aa) 


4ne“mAuR 
uw is the Fermi energy of aluminum and is approximately 
11.6 electron volts. As mentioned previously the area 
common to the two metal films that make up the junction 
1s approximately l mm? . However, as we shall show in 
Chapter VI, tunneling occurs predominantly near the edges 
of the junction so that the effective tunneling area is 


probably of the order of 1% of the crossed strip area. 
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We take A = 10 mee Thus for a resistance of 50 ohms 


we find that P, = 10°”. 
For convenience we choose a square well model 


with transmission probability given by Eqn. (2.18) 


; One ES, 2m 
i.e. P(E.) = r aoe. exp{—2s Pe 5 


Here %6 is the height of the barrier relative to the 


bottom of the conduction band. The exponential factor 
will dominate P(E.) and we expect that tunneling will 
occur predominantly for those values of Ee which mini- 
mize the argument of the exponent. That is, tunneling 
will occur only for those values of E,. close to the 
Fermi energy u. With this in mind we shall disregard 


the preexponential factors which are of order unity and 
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The integral is easily evaluated and gives the following 


result 
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If we take the barrier height to be one electron volt 
above the Fermi level, we can solve Eqn. (2.23)° £or 


y. We obtain e ‘Y x 107/ 


and this is simply the trans- 
mission coefficient for an electron at Fermi energy 
tunneling normal to the barrier. 

Consider now an electron of wave vector k inci- 
dent on the barrier such that 6 is the angle between 
ke end che, normab banpier.. Then the x=-componenc, GF is 


is given in terms of the Fermi wave number kK, by ky = 


ke cos 8. . Thus 


and we can write down the transmission coefficient in 


terms of the angle of incidence 
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If we define 9, as the angle at which the transmission 
es 


coefficient P(6) reaches k its value at 6 = 0, then 
2 
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CrP o - 
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kos e , 
q exp (-y) 


This has the solution 


30 


en a a ne St 
teh (eS.8) Japa Buloe we SW ‘caer 4 sat" pwuX 

~eness af3 yignite ¢2 etdy bas gr: = ie "s asta a 
«Fine is) Fe AGrIsSle As 403 testost 600 not eis . 


i 7 - 


teisisad od od “Esircon paiisa aie 


-louk a 163959V Svsaw to aoxsaste As Wor: tebieec9 | -. _ 
, 


asawied sipns sh at 6 dott dope aeksyBa, ads no 9 vol 


2 26 tastogmep—x eit nodT «ssitszed Learton ors bas A 


= 4 yd .4 asdmun’ svew inset aia 20 uaa iad noviz Be 
é 


@ 


i ! a re . 
ener +8 so ta 
a - 
. § 
. a # 
; BSCS A = Yi 8ob) ae = 
5 ‘ ' . 


7 4 : - = : : : ~ 
a : Ar TAeToOIIse> nobeseimansa3s Shit «awS6b stietw Anus aw bn 
i j - 7 = bh 


} = ; Ne eh: 


o F.* 7 ay i 5 j 
Pert | @uaabioay . sabns od3. io ered 
: — 4 » 7 - : ~~ teal 
ois G yt etn —F 
| | bah per: 14 “a? ‘ge BOS rT] “he F: 
: ; HeOF Fs ze “mH s« (8)4¢ dene 


: 7 i . > 
+ : 
| . colanigensst ads Hodkew +E alias ety Tel a® anitsb td 


me 
- . ‘teas, 0 = a: 46 Sufew atk 2 padese: (874 Jsasiox: ta "I 
- ear Ps a atte. - 


ou 


n 
h- 
+0} 
D 
22 


and fer me l eV, we obtain 6, » 5°. It follows then 


1 
“ 


that tunneling is confined predominantly to those elec- 


trons which have velocity nearly normal to the barrier. 


224 Nonlinear Characteristics of Tunnel Junctions Due 


to Barrier Dependence 


Although the tunneling current is generally found 
to be very nearly ohmic for low bias, this is certainly 
not the case for high bias. A typical example is given 
in Fig. (2.3) where the I-V characteristic of an Al-A1,0,- 
Al junction is displayed. At voltages comparable to the 
barrier height, the transmission probability will become 
a strong function of bias and thus contribute nonlinear 
terms to the tunneling current... Unfortunately even at 
the low voltages that shall concern us in the treatment 
of the zero-bias-anomalies this gross nonlinear background 
is still in evidence. Therefore; we briefly consider the 
effects of the bias dependence of the barrier in order 
that we may separate it from other processes occurring 
in the junction. 

We start with the basic expression given by Eqn. 
(2.13) for the current density. However, as we saw in 


the last section, the current flow across the barrier is 


a 
t bes. 


: : _ - 
revs awollo® o1 ae i ist do “sw ms f AF zor | 
: ea 
-pal<= saor?d at y fr aertimobs3 denddnee ei es a peat 


me © 


_taivted ety oF fewrted yitsean ystooloy iit dio iitw ae 


{ i 


- Sg pabktonul tenn’ Jo epkate ttstosisi seonitaoM § 
- ; sonabysged taizise ea 


a 

a 
a = ¢ . , a 
baitio’ yilstened ef Jnsrivo phileansy sit Hpsaunee 
vieistires ef elit \ en ta wok 10% Simao yitsen ev ® 


zi, eTqaibxs tsotaqyt A agid dgid 162 seso id 


| gi fa-ta & to Sisee ined oasaulp vA at avery " oh) oa 
| | a 


edt ot sihdeteqnop 2er od Lov. 3A inoystigess ai estan sv: 
| impeveck Dietes yo. Ligon ABLE ai wsAy 4s bytes atpiad wobtaz 
= , a I . 
tesrni lion’ -st Jide ies: eons ins BB Ro no rout bali: 


. 36 devs yledonuaaoInl” -Japetbe gaxlost sits oF 


: 
; dremgsexd os ni ev a ears i; ted aspetiov ae 
; J : , 7 { 
; bitoxpased aéacilaon aay 1p jake vj fenone-aassrores ott to 
: Nd a a 
Le sit tsbtenos ylisind ow as otaxkatt” ‘gootishive. «i iittes 
_ _— ’ 


gebio mi torsisd sit to a Sohobasgab este ats to esos ot: = 


: AT | PL 
a ‘pitkaabene (esieaeodsel dol ve : aaa: a aieege yom aa saclt 
tt 7 7 


: en, | «tres ” 

oo aes a vase : - ee | 
< yd. bp 410fses1gKe sid w) ded 
Solcando-enrspeeae maaan! 


aw an | | vam sel xo (CE 


7 oe 


nomen ie 


-suotzoun£ qeuung tw-o°tw-Tw ue ZOJ saazno A-I (€°7) ‘Sta 


3zZ 


Moff =1L (3) 
M008 = (q) 
4oS6C=5 (2) 


0G GE O00v- C0 Ss 


(AW) e6beATOA 


ochre) 


(wu) 


predominantly due to electrons whose velocity is almost 


normal to the barrier. Therefore it is convenient to 


ES: 


use the variable %] = { = ES which is very small compared 


to the Fermi energy. The current density is now given 


Vae2) 


by 
J = anem | P(e)N(e)de (2.24) 
h 
Ul 
where 
1s es 
N(e€) = kT Ln —_ €-eV/k,T . 
l-+-e 


Following the work of Stratton (1962) we use the WKB 


approximation for the transition probability given by 


Home €2019.).4 
X5 2 
i.e. P(e) = exp|-a| (oO(x Vv) = Ww 46) - dx 
1 
where 
x. 2m, % 
Q = Z (=) ° 
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We expand the argument of the exponent in P(e€) ina 
Taylor's series for small ¢€ and retain only the terms 


linear ine. Thus we have 


P(e) = exp[-b(V)-c(v)e] (2.26) 
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where the bias dependent constants b and c are defined 


by 
x5 : 
fom ty Fue) | oxy Aba itle dix |, 
x 
a 
Kee) 
a 72 =e 
eV) = 3 | Lou We =e” dices 
ai 
Eqn. (2.24) for the current density now becomes 
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It is possible to do the above integral in closed form 
with the aid of the substitution 
E/K,T 


ZS 


and the result is 


Teka is 
jesesten sorte ee? aie ee Ce wee (2.29) 
hc sin (mck,T) 


At a given voltage the temperature dependence of the 


current is then given by 
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For 0°K and small bias we note that 
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For a square barrier one can confirm from the definitions 
Of Dang cect vernsby eligi. (2.\cu)eccnat hon. u( 2662) scor SS 
ioeidentical to that ol Bqnaw2.23) where y= b. “Li swe 
use the estimates made of e ’ and P in the previous 
section, we can estimate the magnitude of c and hence 


deduce the temperature dependence of the current 
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Since kp am. 066 meV/°K and ™ yf selL.6 <eV;“the temperature 
factor becomes appreciable only at room temperature and 
can be disregarded at low temperatures. This crude 
estimate is somewhat justified by the I-V curves in 
DELO ts 2hrio)ee 

Stratton (1962) has shown that Eqn. (2.29) may be 


simplified further by expanding b(V) and c(V) in a power 
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series 1n V. He obtains the power series expansion for 


the current 


e 5 
LL = 
pale ae OV Zt dV ee 


where T4 is the zero bias conductance and a ee 
are constants that depend upon the geometry of the 
barrier at zero bias. He further shows that the even 
terms “in this ‘expansion are zero if the bdrrier as 
symmetric and one has then I(V) = -IW). However, for 
asymmetric barriers, he points out that this is no 
longer the case and asymmetry in the conductance curves 
is possible. One finds experimentally that this is the 
rule rather than the exception. Both Rowell (1969) and 
Hartman (1964) have discussed tunneling through asymme- 
tEVO barriers < 

We can summarize this section by saying that the 
ideal current through a junction should be very weakly 
temperature dependent at low temperatures, and it should 


be ohmic at small bias with a small amount of deviation 


present from the non-linear terms. 
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2.5 Perturbation Formulation of Tunneling 

Up to this point the phenomenon of tunneling in 
simple normal metals has been described as the trans- 
mission of electron waves through the oxide barrier. 
In this regard the theory appears to explain the 
experimental details. It was found necessary to take 
a somewhat different approach to explain the observed 
results in superconductive tunneling. This alternative 
approach, advanced by Bardeen (1961), is based on time- 
dependent perturbation theory. The tunneling process 
is regarded as a transition of an electron across the 
oxide with the barrier itself acting as the perturbing 
potential. The electron (or quasi-particle in the case 
of a superconductor) densities of states then enter the 
problem explicitly and this is essential to explain the 
experimental results in the case of superconductive 
tunneling. 

fe merpnot Gitricule te seeathat Chia is oa valid 
description. We first imagine, as in Fig. (2.4a), an 
unperturbed system of two metals separated by a suffi- 
ciently thick oxide layer that no tunneling occurs. 
The total Hamiltonian for the system then contains one 


term, H which describes the left hand metal, anda 


Qe 
term Hs which describes the right hand metal. Each of 
these has a one-electron series of eigenstates. Let us 


denote one of the eigenstates in the left hand metal as 
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Yo and one on the right as es These states are cons- 
tructed to decay exponentially into the oxide. If we 
now decrease the oxide thickness to the point where 
tunneling is possible, (Fig. (2.4b)) then ve and vo 
will overlap. We may now seek a modified solution to 
the time dependent Schrdédinger equation by taking linear 
combinations of be and Wie The transition matrix element 
of the Hamiltonian between states Ve and v, can then be 
evaluated using the usual methods of time-dependent per- 
turbation theory. 

Bardeen (1961) showed, by following this prescrip- 


tion, that the transition matrix. .element M between the 


Yr 
normalized states bo and v is given by 
n? * x 
My. = - 35> | {Wo (x) Vb (x) - WL (x) VY, (x) F.dS . (2.33) 


iS) 


The integral is evaluated over a lamina S parallel to, 
and within the barrier separating the two metals. 
Knowing the matrix element for tunneling between the 

two states, we can immediately write down (using the 
Fermi "golden rule" of first order time-dependent per- 
turbation theory) the transition rate ee for an 
electron in state Po to tunnel into the right hand metal 
fe 


6 (E, - EE). (2.34) 
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Considerpeforvréexampley (the case iofca rectangular 
system where the left hand metal is represented by a 
potential well of width Ly (the thickness of the metal 
film) and the oxide barrier is of constant height Po 
above the bottom of the band. The boundaries of the 
metal are denoted by coordinates Xo Ly and X,) as in 
Figet(2.4a).  A-simple)calculation suchiassiseqiven in 


Schiff (1955) shows that the normalized eigenstates of 


such a system are given by 


one Mee nate Nair pear es eee 
mee AL, sin[k, (x-x,) 1 X,-L)<x< xX, 
hs Gk: aisS -K, (x-x,) 
Ee + is s sail ieaeb 2 Q i mse. 
Ate Ce 3+ ke 1) Q 
bat as Xo 
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where 8 is some constant phase, while the energy of the 


state is given by 


) (23,36) 


relative to the bottom of the band. ky is defined by 


xX 
i, = > 16, — 8, 1, where) BL = EY 427) 
A a Q Q 
The factors in front of the wavefunction are determined 
by the requirements of normalization and the continuity 
of the wavefunction at the metal-oxide interfaces For 


Ly eX 
the parallel part of the wavefunction, e =I ely we have 
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chosen to impose cyclic boundary conditions and nor- 
malization over an area A parallel to the barrier. 
The normalization of the x-part of the wavefunction 
can be restricted to the well since the thickness of 
the metal films is much greater than the penetration 
depth of the wavefunction into the oxide. The wave- 


function can then be characterized by the spectrum 
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where n_, n and n take on integral values. Ina 
SU 29 

completely analogous fashion, we can write down the 
wavefunction vy which decays into the oxide from the 
right hand metal. 


The matrix element between states k, and kee which 


z 
iésjust the overlap integral inwtgn. (2.33), '1s then 
given by 
“2 che cee —K (X)-X)) 
M =- = 6 ees Ooo ie 
k 2m “k k Eppa e Agee Bao 
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Ls eee { noe % xX, 
(23 9) 


The conservation of parallel momentum demonstrated in 
the delta-function is due to the overlap integral in 
Eqn. (2.33), and we have imposed conservation of energy 


so that K a= ee lll If the bottoms of the bands are 
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separated by an energy Eos then 


A-k? 
*r 
E = = EF. + E : (2.40) 
x) 2m fe) Xo 
Since rs ir (the oxide thickness), the square of 


the matrix element may be expressed in terms of the 


transmission coefficient as defined by Eqn. (2.17) 
Dy ee : Ky 
eee ee [7 = eo ogee gee) sat ae. ae 
=U eal ae x Ones 
v 


This may further be expressed in terms of the one dimen- 
sional density of states normal to the barrier (i.e. the 
number of aa states per unit energy range). According 
to Eqn. (2.38) the density of states (which we denote 


by Py ) is given by 


a ee ee (2.42) 
ra 


P(EL ) 
IM. x “= 5 oh y a tee : (248) 


This result, although derived for the case of a square 
barrier, holds in general for arbitrary barriers such 


as the WKB type (Harrison (1961)). 
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We may now proceed to calculate the transition 
rate Posy using Eqn. (2.34). First we must modify 
this equation by inserting the statistical factor 


) fe) ' : 
k (1-f, ) into the summation to ensure that there 
—2 ash g 


ne 
is a State ko available to make a transition into a 


state k which is vacant. Thus 


O O = 
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To this we must add the transition probability of the 


reverse process (i.e. k, may be filled by a state from 


fi 
thesrnight).- Thus 


The total transition rate is then 


ai 2 fe) fe) a 
Sy lo | 2a(fos = for OME ja (2.46) 
Ko n ioe ky k Ls k- ko k 
The summation over ay is done using the momentum con- 
i 


serving delta function in the matrix element while the 


summation over ky is carried out by using the density 
i 


of states and the energy conserving delta function. 
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wi we use Eqn. (2.42) for the density of states and 


remember that the x-component of electron velocity is 


aN hk. /m , then we see that 
= 
Q if 
Ve Raia 5 : “ 
We Vela US es (2.48) 
Ko oo Ke ie 


This result is exactly what we expect since the transi- 
tion probability per unit time for an electron on the 
left ‘to tunnel into the right is just the number of 


times it strikes the barrier per unit time (v, /2L,) 
» 
multiplied by the transmission probability on each 


approach. 


The current across the barrier is just given by 


2eA 
I = 2e ) T = 28! ak [ax (oye é (2.49) 
eS Eee ake) eam yess 


We use Eqn. (2.47) for Ty and remember that E. = 
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E, + nok) /2m CO WET term (in, tie rorm 
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This is the same result for the tunneling current as we 
obtained in the one-electron picture in section 2.1. 


We note the conspicuous absence of the density of states 
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factors in the above result. This comes about simply 
because the square of the matrix element is proportional 
to the particle velocities (Eqn. (2741)] or tnversely 
proportional to the densities of states [Eqn. (2.43)]. 
In the case of super-conductive tunneling it will be 
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CHAPTER] 2 OL 


THE ZERO BIAS ANOMALY 


Ge Nonequilibrium Tunneling 


Wesnoted) im, thes iast chapters (Ban ..(2.45) J athat 
the transition rate for tunneling was given by 
v,P (EL) 


el. Seat ae ace ° 


This can be interpreted as the inverse of the time, Tp 
that it takes an electron to decide to make a transi- 
tion through the oxide. Since tunneling occurs predo- 
minantly at directions normal to the barrier we can use 
as estimates for he the Fermi velocity Ve which is of 
the order of ies cm/sec in typical metals. We saw in 
Section 2.3 that at the Fermi energy yp, P(t) = poets 
Taking the width of the metal film as L x 2000 A, we 


obtain the estimate 


Te = = ~ Ome seconds. Co 


In conventional treatments of the tunneling process 
it is always assumed that the metal electrodes are in 
their equilibrium state. This is accounted for by using 
the equilibrium Fermi-Dirac distribution functions sige 
However, it is obvious that this is not exactly the case. 


When an electron tunnels into the opposite metal it does 
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so by conserving energy and therefore finds itself at 
some energy above the Fermi level in the metal it 
has tunneled into. The electron is thus in an excited 
non-equilibrium state and will rapidly thermalize to 
the Fermi level by a number of possible relaxation 
processes present in the host metal. By using the 
equilibrium distribution functions in tunneling calcu- 
lations it is then implicitly assumed that the excited 
state decays immediately or has a zero relaxation time. 
Nevertheless, for a short time, the electron is in an 
excited non-equilibrium state, the short time being 
characteristic of the relaxation times in the metal. 
While it is in this excited state, it will prevent any 
other electrons from tunneling into that same state by 
virtue of the Pauli exclusion principle, thereby eftec- 
tively “blocking” that state for a short time. 

Consider the situation if the relaxation times in 
the metal were of the order of the tunneling times, 


6 seconds. We would then 


that as, .of the’ order of ite 
expect the electrons that occupy the excited states, 

by virtue of their long life times, to strongly block 
or impede any electrons that want to tunnel into the 
same states. This situation would manifest itself by a 
strong decrease in conductance of the tunnel junction. 
In typical normal metals, relaxation times are of the 

y pala 


order of 10 ~ to 10 seconds at liquid He? temperatures. 
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Therefore we would expect to see deviations in conduc- 
tance due to such blocking effects to be of the order 


of the ratio of the relaxation times to the tunneling 


oe +5 


Limes; that 1s to one pare in 10 10 Small 
changes in conductance such as these could easily be 
seen experimentally by using the derivative measurement 
techniques outlined in Section 1.3 which are capable of 
resolving changes in tunnel junction conductance to one 
part in ee 

To illustrate this idea in more detail, we shall 
briefly consider only the interaction of the electron 
with the lattice. It is well known (Wilkins (1968) ) 
that the relaxation time of an excited electron in an 
electron-phonon gas decreases as the excitation energy 
of the electron above the Fermi level increases. Con- 
sider the situation at 0°K where the Fermi levels are 
perfectly sharp. If our junction is biased by only a 
small amount, the tunneling electrons that are injected 
into the metal have low excitation energies above the 
Fermi surface, therefore their life times are relatively 
long and thus we expect to see the minimum conductance. 
As we increase the bias we increase the energy of the 
injected electrons above the Fermi level thus decreasing 
the effective relaxation time. Since the tunneling times 


remain essentially constant within the bias energy range 


we would then expect to see an increase in conductance. 
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Qualitatively then, an effect of this sort should 
manifest itself as a small dip in conductance of the 
Order sOr Onerpart win ae to lone and a few millivolts 
wide (since this is typically the order of the Debye 
energy). Small conductance Minimavwith this behavior 
near zero bias have been seen experimentally in a number 
of stunnel junctions, by ‘Chen “and Adler (1970)) and others 
(see for example: Rowell, McMillan and Feldmann (1968)). 
In what follows, we shall attempt to Justify this effect 
theoretically and compare it with experiment. 

To treat the blocking effect quantitatively, we 
take the point of view that the distribution functions 
to be wsed in the tunneling current are steady state 
rather than equilibrium. We use a transport model based 
on the Boltzmann equation to find the steady state dis- 
tribution functions and thus their deviations from 
equilibrium. Thesmodel wis wellustrated in Big. (3.2) 
and to avoid cumbersome notation we simply refer to 
Guanti ties on the right: with “ayprime~. 

The rate of change of occupation fun of a state 


Ky on tie right secue {ce athe following processes: 
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is) jUStathe transition rate rgiven by Eqn. (2.44). 
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(2) Electrons can leave the state k' by tunneling 
into the left as in (b) of Fig. (3.1). “Thus 
ofhy 
= 271 2 
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ener Tuer eer Dey 
C3) 
(3) Electrons may be scattered out of k by various mecha- 
nisms. For instance, inelastic scattering with emission 


SleavSOrptionyor, a Lattice phonon vas in (da) of Fig. (3.4) 
will alter both the electron's energy and momentum. 
Elastic scattering by impurities and lattice defects (c) 
as well as non~specular boundary scattering (e) will 
change the electrons momentum. We must remember that 

we have normalized the electrons to a box of thickness 

L so. that specularascatteringtat the sidesiis, implicitly 
assumed. We shall account for these scattering processes 
by introducing a relaxation time Ther! 4 which insgeneral 
will be energy, momentum, and Honoeso une dependent. 


Thus we shall assume that 
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where Atl. is the departure of f 


k! from equilibrium. 


Thatmi $3 


ee 


(3.2) 


where fei is the equilibrium Fermi-Dirac distribution 


function. Thus if we were to turn off the tunneling 


process any out of equilibrium Ath. would decay to zero 
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SEAT Es 
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We shall discuss the relaxation times T and shall jus- 


koe 


tify Eqn. (3.4) in the next section. 


At steady state, Boltzmann's equation simply says 


that the net rate of change of fi is zero for any value 
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If we substitute the relevant quantities into Eqn. (3.6) 
we obtain the equation 
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into state k' is equal to the rate at which k' relaxes. 
A similar equation may be immediately written, by sym- 


metry, for a state k in the left hand side. That is; 

=f) | lati oe eels @, iol B an (3.8 
FA K! Mkk! k! k k Kk"? . ) 

which just expresses the particle-hole symmetry of the 


problem. Eqns. (3.7) and (3.8) may be reduced by carry- 


ing out the summations as in Eqn. (2.47). We obtain 
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sane el Bape | Mal (ay ee OE) 
Ty sal k k k' 7 
(3.9) 
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P(E.) 
fc ee x (3.10) 
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Making the substitutions a = Ey + ad and = a EeetAtiny 
we are able to solve Eqns. (3.9) algebraically to,obtain 


the difference of steady state distribution functions, 


i, 2 8 , which we shall use to determine the tunnel- 
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ing current. The result is given by 


ia : 3 elt) 


Zi 2 
1+ = 1M PIT 


{Pp Te 


x a a} ye — 
_ - * - vo a 

— | if . a a we 
oe re ia Stee 9 eid of Ladpe ck 5; Aas 
meme yeh God Sonw ona rbone od yysit rus ssup 


J tai tet -vbie bran st sas ot 2 a + . 


. 
Lis 


A (fc) ‘dag - = ga a4 3 -ry) haga rgd 
_ 


7 ‘ 
ait to yxsounva olott-oloizvis 9 ails eseastges sent bes il 

| 
-yits> vd Bessybor sd vem (8,£) fap (3.5) .anpe mete 4 


Eo 
Ristdo oW .(TReS) fp ms 25 ani civamnea ort 40 | pat 


rth 
{ 7 a 7 uch “ re = Paes 70 
og x 4 - at 
(@0E1) 
7A 
" va taht Sey teal. a a 
; a — at x a 
; Ln 
(sa) ig te 
pee aa a ee 
" a ath 
7 | | a) - Ww i ri : 7 
7 4 F taal 1 S143 pris Pen as ; Pe: = PE anottust 2 ert ie 


ataldo ot yite irweauith: (ey aes aries 2 otds ois, § aa 


vaio: Pail riick sik a2 id i sive fy stripes ts sudiersI tthe 


- 


stead oa ET sian = = , 2 2 “ 
7 - a 7 wi nena rnc =. dees 0 pa 
ce 4 

a x Le ee 


54 


Consider the term 


27 Z 
= iM) ey. 
appearing in the denominator of Eqn. (3.11). Using 


Gheyde tins tion iiweniby eqnat use L0)) and\expressangi iene 
density of states in terms of the x-component of velo- 


city, we have 


P(E.) VP (EL) 


Zi [m|? O = 1 a 
Si ep Wantneol ae S fills! 
Al 
= TH (es, 282) 
B 


v . . . . 
where Th is just the tunneling time for an electron in 


the right hand metal. Eqn. (3.11) may now be rewritten 


as 
O fe) 
as 
ne = Se < SL) 
c 4 ae 
1+—-+— 
ve Tp 
Since the tunneling current is proportional to fy. = fn 
we Can see from Eqn. (3.13) that the deviation from the 
ideal Current Will be of the order of T/T. , the ratro 


of relaxation times in the metal to tunneling times. 
We now turn to a discussion of the relaxation times 


and their behavior with temperature and energy. 
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3.2 Relaxation Times 


We shall describe the scattering due to lattice 
defects, impurities or non-specular boundary scatter- 
ing by a constant relaxation time T;- We suspect these 
processes are elastic in nature and should therefore be 
temperature and energy independent. Since we lack | 
detailed data about the structure or impurity composi- 
tion of the evaporated metallic films that make up the 
tunnel junctions, we make no a priori estimates of cS 
except to say that it is probably comparable to the 
electron-phonon relaxation time in order of magnitude. 
We shall therefore treat T, as a constant to be extrac- 
ted, if possible, from our model. 

In contrast, we expect that the electron-phonon 
relaxation times will be sensitive functions of both 
excitation energy and temperature. The reasons are 
Stacistical in natures -Atb“0°K,- an ancrease in the 
energy of an excited state means an increase in the 
number of empty states that it may make a transition 
into. ee ltheretore tne probability Olma toast. tion 
increases while the relaxation time decreases. An 
increase in temperature will smear the Fermi level and 
increase the number of states. Therefore we also expect 
the relaxation time to decrease with increasing tempera- 


ture. The phonons, obeying bose statistics can be 
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created in any number. A detailed calculation of the 
electron-phonon relaxation time is therefore in order 
Since by changing the junction bias, we are effectively 
changing the excitation energy of the tunneling elec- 
trons that are injected into the metal. The tempera- 
ture dependence is also important because at a few 
degrees Kelvin the thermal energy k,T becomes compara- 
ble to the energy width of the zero bias anomalies. 

We shall calculate the relaxation time of an 
excited electron in the presence of a lattice (phonon 
gas) by using Fermi's golden rule and considering only 
first order scattering processes in the square of the 
electron-phonon matrix element. In Fig. (3.2) we con- 
Sider the four processes that alter the occupancy of 
Cie exci Led citatenp. Ey is the energy of the excited 
state with momentum hp while Ey is the energy of the 
final state with momentum hk. “The phonon has momentum 
+hq (where gq = p-k) and energy fiw(q,A) where A is the 
polarization index. The processes in the ‘first column 
lead to a decrease in occupancy of p by absorption or 
emission of a phonon while those in the second column 


pb! the 


occupancy of the state p , can be written down by ins- 


lead to a decrease. The rate of change of f 


pection of the diagrams and by noting that an in. baal 


electron or phonon has the occupancy factors f and N, 
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Processes that lead to a decay of an excited 
electron state. 
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respectively, while the final state factors are (1-f) 
and (1+N) respectively. N is the Bose statistical 
factor for phonons. As an example the first process 


in Fig. (3.2a) contributes the following term 


ee 45 2 is enw 
se ae e-pqettr ald fa Elida) (sy fw (q,A)) 


(3.14) 


The electron-phonon matrix element, G(q,i), has the form 


(see Bardeen (1937) ) 


eae ti - 

Slo lees i Soqoon Ly ae 
where MN is the mass density of the lattice and gy is 
the phonon polarization vector. v(q) is the Fourier 
transform of the ion pseudopotential. When all four 


processes in Fig. (3.2) are taken into account we have 
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P 27 2 
ie ae Wee eae {6 (E,- Sag ee Ve, x 
x [£,, (1-£,,) (14Nq) - EN, (158) + (3.16) 
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previous section we let the departure of the state p 
from equilibrium be Ar, , where Af = = £°)\ owe" shall 


use the equilibrium distribution functions for ft. and 


a Thus we are assuming that the electron that has 


tunneled into an excited state p is scattered by the 


Fattree into the ‘equilibrium distribution.” With these 


considerations Eqn. (3.16) now becomes 
ig ie 
= i - 3.27) 
a eee ° 
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where the relaxation time We is given by 
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This expression may be simplified by defining the exci- 


tation energy Ep = ah -—y above the Fermi level yp. The 


summation over k states becomes 
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the energy surface defined by E. = A7k*/2m, and v(k) 


is the velocity defined by v(k) = |V,E,|/f = A|k|/m. 
Since w(q,A) is of the order of the Debye frequency Wy 4 


the dominant contribution in the ¢« integration comes 


k 


from Ee, < Aw, << wu. Hence we can extend the lower limit 
in Eent (3.19) to infinity and evaluate the surface in- 
tegral on the Fermi surface. If we make these changes 
and Canty out the Ex integration over the energy conserv-_ 


ing delta function, we obtain 


Palin Bie sega aoe Vie | sess lGlq ale x 
ee ee ip (Qn) 3 fiv,, (k) x a 


x (2N° (w(qyA) )+14£° (e Hiw (q,))-£° (e,—fiw (g/d) ) } 


(3209) 


We can separate out the phonon contribution in a con-- 


venient manner by inserting a delta function in (3.20). 
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where by definition 


5 ois ¥1G(q,r) | 76 (iw—Fw (q,A)) 
a. (hw) F (fw = YL Gia pr 6 (Aw-fAw (q,A 
= Ee Ci ee ae a 


(252.2) 


- 


i; — re big, eae bins fi ya besitod ove tities, 
a Tutt hae Nag HY micas ya berth an. ep 
ZZ +e yoasups al ayaed ort a6 aebio ait ta ae 


; ; asimos aoitatpetaz # and mi naLtudizdao2 
i ‘okmeL tewol siff Baetee tec =u socom .y >> of spa 


es -ni sostvue oH3 steclave bas ysintiat od (Of 28) ‘ae 


aay 
zopnsis svetd S#amn ew 11 «sortase ipsot edt ae Ieipged 


- « 
~Vi8R Hic? YPIENS sie devo aobterpssni # Stic to YETs - 


7 . 
. 


niatiio ae , aoitonut Bi feb” ord 


c (A), ie .* V as a m : 
7 * “1 {k, 5) > <> =o eo 
aes ar Dage \ Eke, ne ace ae 


7? . PLGA vg) ite 59) “3 -( (4 1p) wits) pete tt <+B) wy CMS) t 
; A, . 2 
ye) ; es 


a : a : } 
~— | 


Ants & nm yotsudti 32 TOD aotewle 3a3 710 | 9PREBYSS, ago OW 


: ~AUS.2) nk notisayt pt tab 5 pe ivagens vd osnnsm ‘seenen 
: . be 
~ ; ; od / 7 > i 
x! ide fem ; Ai(iift) & | fe 9.i 
s = 
7s 7 7 —@ 6 7 7 ae a ‘= Mp - 


2 - as ' 


: 68-8), Lit at tu i 4 (cai Koos & 


7 
od 1? _ 


a: | noite S190 


61 


We are tunneling into an evaporated film which is 
probably polycrystalline and characteristic of an 
isotropic material. Therefore, we do not expect to 
see lattice anisotropy in the relaxation time and 


consequently we average the inverse life time over the 


Fermi surface. .Eqn. (3.21) becomes 
= - = | d (hu) a” (hw) F (hw) [2N° (hw) +14£° (e+fiw) -£° (e-fiw) ] 
0 
(225) 
where 
dS_,(p) ;ds,, (k) 
were —Y__— ¥ |G (q,A) (76 w-fiw (q,2)) 
a eae. Aba wou Me Bian! Yao abictan ents te = 
Qa WwW (OO) 
ae 
Vp (p) 


(Zea) 


The quantity 07 F is dimensionless and while difficult to 
calculate from first principles, can be extracted from 
tunneling experiments in superconductors and is presently 
known for Pb, Hg, Sn and In (see McMillan and Rowell 
(1969)). The awkward notation a” (tw) F (fiw) was invented 
with the following intent: the phonon density of states 


in a metal is given by 


F (hw). =") S 5 | dq 6(fiw - fiw(g,aA)) . (Su25) 
eG) 


Thus a” (Aw) is taken to be a measure of the strength of 


the electron-phonon coupling. For convenience we replace 
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fw by w and refer to w as the phonon energy (typically 
of the order of meV). We shall now proceed to discuss 


the properties of the relaxation time t(e) given by 


5 =r | dw oh enn eR ee S (edu) Se en (3-26) 
0 
where N° (w) = | and Seve epee! Bined (eno 
W k,T e/k,T : s 
e - 1 lt+e 


We first note that the relaxation time at any tempera- 
ture is an even function of excitation energy e«. This 
follows directly from the above definitions. At T=0, 

w) = 0, and the statistical factor in the integrand 


of Eqn. (3.26) iS just) a ‘step. funetion 


fe ft (ee) = Ee) =) =b o<0y <6 
0) otherwise. 
Therefore at T = 0O, 
if Z ie eS 
ee ee ; B28 
=e) = | dw a F(w) , ( ) 
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where the absolute value on the upper limit ensures 
that the relaxation time is both positive and even. 
At finite temperatures the cut-off in the step function 


is smeared by the thermal energy k,T but iS. Still tearay 
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sharp if k,T <<e€. Physically this just means that an 
electron of excitation energy € cannot decay to equili- 
brium by emitting a phonon of energy w greater than ¢. 
Since the excitation energies will be of the order of 
or less than the bias energy we can therefore disregard 
a7 F (w) EOL a> GN ues The maximum voltage we deal with 
in the zero-bias anomalies is ~ 5 mV, therefore we can 
disregard any information in oF (w) beyond w ~ 5 meV. 
In order to facilitate calculations, we shall make use 
of this fact by assuming that the metals we deal with may 


be represented, up to w ~ 5 meV, by a power law behavior 


of the form 
BD < 
Os. Cie Gl) ae te (3.29) 


The relaxation time then becomes 


Eqn. (3.30) may be simplified somewhat to yield 
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and I and ¢ are the Gamma and Riemann Zeta functions 
respectively. Rhodes (1950) has studied the functions 
F {2) and has tabulated their properties and values 
for integral values of n. We shall consider only the 
cases h- = ] and 2. For n:= 1, tle) can be evavated in 


closed form to give 


Ta 
is 7 1 Zz 2 
ae ynenap HP che Ebley ¢ Bee) 
The relaxation time for the case n = 2 cannot be evalu- 


ated in close form. However, Rhodes gives convenient 


expansions for F(Z) with which the life-time may be 


approximately evaluated. For the case n = 2, we have 
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i 27a, 3 
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The Debye frequency spectrum F(w) defined in Eqn. 
(3.25) varies quadratically with w at small frequencies. 
The model with n = 2 then simply states that the elec- 
tron-phonon coupling, a”) is constant at low frequencies. 
With this behavior, the relaxation time should, according 
to Eqns. (3.37),.vary inversely as the temperature and 
energy cubed. Measurements of the relaxation frequency 
in lead and indium by cyclotron resonance (see Goy and 
Gastaing (1972)) indicate that this cubic behavior is 
observed. However, theoretical one-OPW calculations by 
Allen and Cohen (1970) and Carbotte and Dynes (1968) 
indicate the possibility that oF (w) behaves linearly 
with w at low frequencies due to Umklapp processes. We 
shall not attempt to resolve this question here but shall 


calculate the zero-bias anomaly with both models. 


3.3 Calculation of the Current, Conductance and Deriva- 


tive of Conductance Due to the Blocking Effect 


In order to calculate the tunneling current, we 


shall assume that the only relaxation processes of 
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importance are the ones described in the previous sec- 
tion. We characterize the relaxation times for elastic 
scattering and inelastic electron-phonon collisions by 
LA (constant) and Tap (energy and temperature dependent) 
respectively. Since the rates of change of occupation 
of a state due to these two processes are additive, the 
total or combined relaxation time tT is given by 
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for any state ek. 

The tunneling current is simply given by summing 
over all transition rates from left to right (see Eqn. 


(2.46) ) 


KF} 6 (EE (3239) 
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where fe and fis are now the steady state distribution 


functions on left and right sides respectively. Their 


difference is given by Eqn. (3.11) 


(3.40) 


where ae and Tp, 1 are the total relaxation times on the 
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the change from equilibrium to be small, we expand 


Eqn. (3.40) to obtain 


fo) 27 2 
fy Ey = (£1) {1 Se el (Oy. T 1 + Py ty) t- (3.41) 


The tunneling current may now be represented by 
TBS Eb htast (3.42) 


where 
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K! 1) 


is the conventional tunneling current calculated with 
equilibrium distribution on both sides and discussed in 
Sections 2.1 and 2.4. I is the small negative correc- 
tion to the ttnneling current due to the blocking effect 
and is given by 
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where 


We must also remember that the relaxation times depend 


only upon the excitation energies, i.e. Ty a TE oy) 


and Uae T(Ey yo u'). Since the Fermi levels are sepa- 


rated by a bias eV (i.e. u-y'=eV), the expression for 


the blocking current reduces to 
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21ho' (E) 2nfip (E,) : 


In contrast to the conventional current I, the 
blocking current 6I explicitly depends upon the densi- 
ties of state in both left and right hand metals. 
Clearly, if the relaxation times become very small, 
the blocking current $I vanishes as we expect. We may 
note that the blocking current depends upon the square 
of the transmission coefficient and should therefore be 
a sensitive function of oxide barrier thickness. Eqn. 
(3.46) may be simplified by making use of the fact that 
tunneling occurs predominantly with velocities normal to 
the barrier. Therefore we replace the quantity P(E,)/ 


2mfip (EL) by its value at the Fermi energy 
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P(E, ) Py y 


ia aC XS) SC we 


(37047) 
where VF is the Fermi velocity and Tp is the tunneling 
time for an electron as defined in Section 3.1. 


Therefore Eqn. (3.46) becomes 
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is the ideal zero-bias conductance defined in Eqns. (2.20) 
and (2.21). Eqn. (3.48) may be written in a more com- 
pact form by making use of the variable « = E=-y to obtain 
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We can simplify the above equation still more by noting 
that the energy dependence of the effective life times 
comes from their electron-phonon parts which are even 
functions of their arguments. Using this fact and the 


identity 
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the blocking current may be written as 


io} - : 
st (vy =e | de [£° (c) -£° (e+ ev) ] (Ee) + te). 
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The conductance of a junction, normalized to zero 


bias is given by 


Cl eo ee) 
haaie a => Gm lau AF ; ee (33,4) 
fe) fe) e) 
where 60(V) = “ (6I(V)) is the blocking conductance. 


Using Bqn. (3.50) ~we find 
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This is the basic form that we use to compute the cha- 
racteristics of the zero-bias anomaly. It may easily 

be verified from Eqn. (3.52) that the blocking conduc- 
tance §60(V) is an even function of V and hence it is 
symmetric about zero bias. In subsequent work we shall 
simply use V in place of eV and define the bias in units 


of energy (typically meV). 


Zero Temperature 
In order to observe the behavior of the blocking 


conductance §o0(V), we set T = 0 since in this case we 


Soneioybhos paitsold srt ei ((V)15) os = (Vj vt sxedw 
bait ow (02.6) apa pakeal 


eee, i ES, fad ae van |= 


Na aii Maltigucat adeet aa ee na 
yliebs ysat +1 OLE, AEE Oe ES 

-cuhdas enitoblil/ads InAs : ($2.6) Lop nosh bal apr ad 
ak dE sonal bas Vi fo noissay? nove ms et ins 


Fig. 


(os) 4 


Conductance and derivative conductance 
for n=). 


eV 


oe 


ea Tawa 


Pig. (3.4). Conductance and derivative conductance for m2. 


73 


can evaluate it in closed form. At zero temperature, 
fe) 
-df- (e)/de = §(0), and Eqn. (3.52) becomes 
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where for convenience we have assumed a symmetric 


tunneligunction sor that coves «(Vv )ipeand ie = T 


The effective relaxation time 1(V) is given by 
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where Ts and Tep are the impurity and electron-phonon 
relaxation times respectively. At T= 0, Tap is given 
DYE Oies (C3020) 5 
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where n corresponds to the power law behavior of 
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The features of the conductance and its derivative are 


plotted in Fig. (3.3) for n=l and in Fig. (3.4) for n=2. 
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Shawl! Comparison OF TnReCOLY wath Experiment 


In Fig. (3.5), we show the experimental conduc- 
tance of an Al-Al junction at biases up to 5 meV for 
several different temperatures. The dip in conductance 
about zero voltage is what we commonly refer to as the 
zero-bias anomaly. The sharpness of the anomaly broa- 
dens considerably with increases in temperature. 
However, we notice that the conductance over the 
range of +5 meV is not symmetric about the region. 

We assume that this asymmetry arises from the ideal 
background conductivity \discussed in Section 2.4. 


Therefore we suppose that the experimental conductance 


in Fig. (3.5) may be represented by an expression of the 
form 
2 
} = 
or CVs Ty + 04V + O5V + 60(V,T) (3358) 


where the first three terms constitute the ideal back- 
ground conductance (see Rowell, McMillan and Feldmann 


(1968)) with o, and O54 being small constants. One possi- 


i 
ble way to test this conjecture 1s to use the) fact that 
the ZBA conductance given by 6é0(V,T) is an even function 


of bias. Therefore if we take the odd part of Eqn. (3.58) 


we obtain simply 
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and so theoretically the odd conductance is directly 
proportional to bias. The experimental odd conductance 
obtained by using the data in Fag. (3.5) sts sshown an 
Fig. (3.6). We note that it is proportional to bias 
Up cOes meV Co vabout one part in rO7 thus agreeing with 
theory. 

In principle then, we can effectively compare 
theory to experiment by taking the even part of the 
theoretical conductance 


v2 
i.e. Oo, (x) 0s + O5V + 60 (VG L) (3% 60) 


and comparing it to the even part of the experimental 
conductance. Unfortunately we do not have sufficient 
data to isolate the background part ay Certainly 

it should be smaller than the linear term 0,V and we 
shall therefore neglect it. However it may be trouble- 
some in the sense that it masks the asymptotic behavior 
OL the 2ZBA conductance. 

In practice, 4 more Sensitive and accurate measure 
of the ZBA conductance is given by its experimentally 
determined derivative as outlined in Section 1.3. 
Therefore in the rest of this chapter we prefer to test 


the theory by comparing the theoretical derivative of 


even conductance 


; do (V) 
i.e. S = 2 go(v,T) (2.52) 
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with its experimental counterpart. We assume at the 
outset that the two electrodes making up the Al-Al 
junction are similar in nature and have the same cay 
Top and Tye Therefore the even derivative of conduc- 
tivity is given by 


sik do. WV) = A (do (V)) irks eye gf (e)) we ee 
Se, dv dv ome J€ dv Tp 


(3.62) 


For finite temperatures the integral in the above expres- 
sion must be done numerically. The reason is simply 
that the half-width of the distribution =f de 0/ oe is 
About 5.) 5 kT so that at 2°K, the thermal smearing is 

of the order of half a meV. This is comparable to the 
width of the ZBA conductance as may be seen in Fig. (3.5). 
For numerical calculations it is convenient to cast 


EGh.e (3.05) Into ene form 


do (V) 
= 5 = | an —— {S(V-xk,T)+S (V+xk,T)} (3.63) 
70 Ms 0 (ie). (ae as} 
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where the electron-phonon relaxation time LY BS gene is 
given by Eqn. (3.33) for the case where n=1 and by 
Eqn. (3.34) for the case where n= 2. 

Fig. (3.7) shows the experimental do ,/dv vs V 
curves for the same Al1-Al junction measured at T= 3.0, 
4.1, 5.0 and 6.0°K. The dashed line is the theoretical 
curve given by Eqr. (3.63) for the case where n=1 (i.e. 
o7F (w) = a,w). The theoretical curve was fitted to the 


3°K experimental data and yields the values 


ee =e dicots eens 


285 <0 4 seen 


a 
I 


The electron-phonon coupling constant a, was chosen to 
agree with the calculation of oF (w) F (w) by Carbotte and 
Dynes (1968) and was given the value a,= -003/meV. The 
overall shape of the theoretical curves is not extremely 
sensitive to temperature. Accordingly, the theoretical 
and experimental peak positions and heights are compared 
in Fige (3.8) forg the Same setyof values Tar Tp and aj 
Excellent agreement is found for the peak positions as 

a function of temperature. The theory predicts a stron- 
ger temperature dependence for the peak heights than 
experimentally observed. The behavior of the theoretical 


curves are found to be quite insensitive to the choice 


OF Ts. Fig. (3.9) depicts the temperature dependence of 
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Fig. (3.9). Dependence of peak position on Tae 
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the peak amplitudes with different T'S. We have fixed 
the theoretical amplitude to agree with experiment at 
3°K. Fig. (3.10) shows how the peak position behaves 
with changing Tj. We note that there is relatively 
little change in the characteristics when the impurity 
time T, is varied over an order of magnitude. We may 
conclude that the zero-bias anomaly shape depends weakly 
upon the nature of the elastic scattering mechanisms in 
che. faim. 

We have also calculated the behavior of the zero- 
bias anomaly for the case where n=2 (i.e. ae aye = 
a5 aie In Fig. (3.11) we show the second derivative of 
current for n=2 at 3°K. For comparison we also show 
the experimental data as well as the n=1 result that 
was illustrated in Fig. (3.7). The n=2 model curve was 


fitted to the experimental data as shown and yielded the 


values 
23) 
aes Spee 40) sec 
7 x10” se6c 
Tp = ¢- 
where the electron-phonon coupling constant a5= s00037 


mev- was again chosen to agree with the calculation of 
Carbotte and Dynes (1968). The theoretical and experi- 
mental peak heights and positions of the second deriva- 


tive are compared in Fig. (3.12) and PQs Sco) 
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respectively for different choices of Tse The results 
again show that the peak heights and positions are 
relatively insensitive to T,- In general it may be 
noted that the overall behavior of the n=2 model is 
very Similar to that of: the n=]\ model, although the 
agreement with experiment is poorer in the case of the 
former. 

An additional example of a tunnel junction that 
displays the ZBA behavior is that of a Pb-Pb0-Au junc- 
tion whose experimental derivative of even conductance 
is shown in Fig. (3.14). Although’ the ZBAwstructiires is 
Clearly in evidence at low bias, it is masked to some 
@xtent by a rising background. This is probabisaidue to 
inelastic tunneling of electrons at the lead-oxide in- 
terface. In order to account for thisein a theoretical 
calculation, we have simply fitted a linear background 
to obtain agreement with the experimental data at 1.13°K. 
The theoretical curves at higher temperatures were then 
generated with the same set of parameters. The calcula- 
tion was done using an n=l model. We find that although 
the shapes of the curves agree at higher temperatures, 
theory predicts a stronger temperature dependence of the 
amplitude than is measured experimentally. 

In conclusion, the observed structure near zero 


bias in the conductance and its derivative do/dv for 
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normal metal-insulator-metal tunneling junctions has 

been accounted for by including the effects of finite 
electron relaxation rates in the junction electrodes. 
The non-zero life-times lead to an observable block- 

ing of otherwise available electron tunneling states 

at low bias. The fit to our data suggests that 

a” (w) F (w) behaves nearly linearly at low energies 


rather than quadratically. 
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CHAPTER IV 


INELASTIC ELECTRON TUNNELING 


4.1 The Inelastic Current 


In our treatment of the zero-bias anomaly we 
assumed the oxide barrier is without structure and 
provides essentially a homogeneous forbidden region 
to the electrons as would a vacuum. However tunneling 
measurements have shown that a considerable amount of 
non-ideal behavior is present in I-V characteristics 
which may be directly attributed to the presence of 
the barrier. A typical example is the Pb-Pb0-Pb junc- 
tion whose experimental second derivative of current is 
shown in Fig. (4.1) as a function of bias voltage. We 
note the large amount of structure present in the second 
derivative up to 60 meV. This behavior is due to the 
fact that electrons may lose energy as they tunnel 
through thefbarrier. The oxide is host to a number jor 
mechanisms that may cause the electron to tunnel inelas- 
tically. The oxide may have bulk phonon structure as 
well as ever present impurities that provide molecular 
rotational and vibrational excited states. The metal- 


oxide interface may also provide additional degrees of 


freedom. 
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If a voltage V is applied to a junction, the 
tunneling is elastic to a first approximation and the 
I-V relationship is linear. However, suppose that at 
T= 0 there is a system in its ground state that re- 
quires an energy hw, to be excited. If the bias energy 
eV is less than fwys the electron cannot tunnel inelas- 
tically by exciting this system because there are no 
available states in the other metal. When ev = Aw, ’ 
the electron may now tunnel inelastically. Thus an 
electron in a definite energy state E has more than one 
available energy state into which it may make a transi- 
tion. The opening of this new threshold will result in 
an additional current and therefore a small change of 
slope in the I-V curve at eV = fw: Consequently there 
will be a step in the conductance and a spike in the 
derivative of conductance at this bias. In practice a 
finite temperature will smear the Fermi surface of both 
metals and will thus smear the spike although it will 
be centered about eV = hw: In general a continuous 
excitation spectrum will be present in the oxide. The 
second derivative of current will then mime the frequency 
Spectrum of the barrier. In principle) theptunhedegunc- 
tion will then act as a spectroscopic probe of the 
barrier and any impurities within it. Lambe and 
Jacklevic (1968) confirmed in a series of well prepared 


experiments that this was true in practice. They showed 


yeishe eei7 ots Ti {Sad tors od ae 
-ealani Lennny Joanse, HOTELS ont gut ate ore 
on Sts sedis sated mageve aids pee ya te 
\ gut =) Ve rete fst ont sadso ailt wi estnza ottel dave 
16. auat -yilsoiseslont iénauts won van iettebte aria 
ano asa} siom es 4 stci2 ypsseas estinitab 6 at alata 
-tens3x3 2 Sia yem di noid oshl, state yoisas otdeziovs 
“i aivess Lfiw Biodesri13 woh aiid to patrhsge sat ‘aoks | 
to spnéAo ‘fteme 5 Siotexsiy bas josaigo isnotsibbs ne 
stent yisasnpesnoo «wit = Vs 38 S¥auUn V=-I sd at aqole 
ait ol oitqe 6 Bas Sonstoilbneso ont ai gate a sd Iliw 
& 9oij201q ai .asid eidd Is sdristouhacd to ovisavizeb 
‘Zod 36 sostive imasa sda xrseme Iliw sruteséqmiet sdinit 
ifiw +¢ dovodals tee! sey) teenie code Dbw bas efsjem | - 
auounitnos s [atesnap nt ..wih = Ve gueds Bersinss ed. 
vit ,Hijx0 sis ai Jnbssig r 1 Sima ozsasta 
erent od omim osit Liitw Secane to mrigavieed baocose 
eomdt Larnst sd aks" — sft Io muxjoeqe 
ort 40 scoxq vigooe 
anaes ha 


94 


that one could observe and identify molecular vibrations 
in doped M-I-M junctions. They also showed how one 
could estimate the effect of temperature on the line- 
width of a spectral peak due to the smearing of the 
Fermi surfaces in the metals. Weesharl obtain their 
resultsherecin a slightly different fashions 

The mechanism whereby coupling occurs between the 
tunneling electrons and internal states of the barrier 
may be thought of as a perturbation of the barrier po- 
tential. Consider the transmission coefficient for the 
WKB barrier given by Eqn. (2.19). If é6(x) is some small 


perturbation to the barrier potential 9(x) we have 


x 
re 1. 

P (B,.) = exp|- | 3 [> (x)-6(x)-E,]}* dx . 
ae 


If 6(x) is small we can make an expansion to first order 
im ox) to. obtain 
Xx x 


Ae 2 2 2 
2m ” Pe aes 
P(E, )= exp(-| bale )-EL J] ax [2+3 | ge b(x)-5,) ax 


it) 5 § 


The first term will give rise to the elastic current 
while the second term will give the inelastic contribu- 
tion. However,this is a very over-simplified view of 
the inelastic process. The perturbation potential will 


be the result of a superposition of potentials which 
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depend upon the positions of the oxide atoms, the impu- 
rity molecules and perhaps the positions of the metallic 
atoms diffused into the oxide. As these positions are 

a function of time, they will contain factors like 


SMCS is) where Wg is the frequency of the normal modes 


of wave vector q of the system. Thus if one were to do 
the problem in the context of a time dependent pertur- 
bation expansion one would obtain energy and momentum 


conservation equations of the form 


oes: = le ; DEO DK" = fg 


whereyEy Vand hk are: thesenergy andgmomentumforethe elec- 


[~ 


tron in the left hand metal and E and Hie. are the 


k! 
energy and momentum of the seen abter=1c-has tins 
neled through the barrier into the right hand metal. 
With this in mind we can treat the process of inelastic 
tunneling in exactly the same way as we determined the 
electron-phonon relaxation time in Section 3.2. That 
is, we suppose there exists, a matrix element Mt that 
connects the states k and k' on either side of ane 
barrier inelastically. If we make use of Fig. (3.2) 
and imagine that a barrier separates the initial and 
final states k and k' respectively then we can write 


the transition rate for the process whereby an electron 


in state k on the left tunnels inelastically into the 
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right hand metal. That is, 


eT 2 fe) 
ES LIM | (6 (BB y Hw) (Ey (lf 41) (L4NQ) “£9 1 (1-8, NG] 


7 fo) OV OLD fe) 
+ ae Mae OP ee rene Ng ei a) 


(4.1) 


where a7 and fe 


the left and right hand metals respectively. NG is the 


are the equilibrium distributions in 


equilibrium distribution function for the internal degrees 
of freedom in the barrier, and we assume that Ng depends 


only on ee We assume that the left hand side is at bias 


energy eV above the right hand side. Therefore 
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lte-r- l1l+en- 
We define oe = fpf i (V=0) and make use of the fact that 


the transition rate is zero for zero bias. Therefore 


Eqn. (4.1) becomes 
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according to our density of states and box normaliza- 


tion scheme presented in Section 2.5. Carrying out 


the summation in Eqn. (4.3) we find 
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where we have defined 
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The inelastic tunneling current is simply given by 
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ne a 
él = | de [£° (ce) -£° (e+ev) ] | dw G(w) [2N° (w)+£° (e+w) 
—% 0 
+£°(w-e) ] , (4.9) 
where 


We have expressed the current in terms of the dimen- 
Sionless function G(w) defined in Eqn. (4.8). We note 
that G(w) is defined in a fashion analogous to the 
function ow (s) as in Eqn. (3.24). As such G(w) may 
be thought of as the product of the spectral density 
of states of the barrier multiplied by an effective 
frequency dependent coupling constant. It is not our 
intent here to attempt a calculation of G(w) from 
Eqn. (4.8). However it will essentially be dependent 
upon the transmission coefficient and the dynamics of 
the coupling of electrons to internal degrees of freedom 
OL the barrier. 

We shall now assume that N- {w) is the equilibrium 
Bose factor as defined in Eqn. (3.27). Thus Eqn. (4.9) 
would describe the inelastic tunneling current due to 
the bulk oxide phonons in the barrier. We note that 
because of the properties of the statistical factors 
Uieeei 49) Ene inelastic tunneling current is an odd 


function of bias voltage. Thus the conductance due to 
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it will be an even function of bias. This has been 
established experimentally in the case of Pb-Pb0-Pb 
junctions by Rowell. McMillan and Feldmann (1969) 
and is true in general for any junction where the 
oxide phonons may be identified easily. We also 
note that the T = 0 limit of the inelastic tunneling 


current as given by Eqn. (4.9) becomes 
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Therefore the inelastic zero temperature conductivity 


is given by 
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and thus 
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For finite temperatures, Eqn. (4.9) may be sim- 


plified considerably ‘by "reversing the, order of eintegra- 
tion. The integration over the statistical factors 


may then be carried out analytically and yields 
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a5 (62)"="c | dw G(w). a KB) - KR (4.14) 
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where the kernel function K(x) is defined by 
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(4.15) 


ands Ss an“even: functionsof; sesin weqn jo 4s0h4)e thesterm 
K ((V-w) /k,T) represents inelastic current flow by 
emission of a phonon of energy fiw whereas the term 

K ((V+w) /k,T) refers to inelastic current flow with 
absorption of a phonon of energy fiw (anti-Stokes flow). 
If we have an extremely narrow spectral line such that 
it may be approximated by a delta function distribution 
(apse. Gili) d(w-wo)), then the second derivative of 


current due to this resonance is given by 


O 
S (61) = Ss Kee )- KG) (4.16) 


ih thespositivonser: the mesonance;isymuch: greater, than 
the thermal energy (i.e. Aw, >> k,T) the absorption or 
anti-Stokes term is negligible and we are left with 
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Using Eqn. (4.15) we find that the response in the 


second derivative takes the form of a dustraibutzon of 
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height proportional to c/6k,T and linewidth 5.44 k,T 
centered at ev = Aw: This is precisely the behavior 
that Lambe and Jacklevic (1968) found experimentally 
and justified theoretically. This analysis describes 
the linewidth expected for a sharp resonance due to 
the thermal smearing of the Fermi levels in both normal 
metal electrodes. Although there may be other sources 
of linewidth broadening in the actual phonon spectra of 
the barrier or the molecular vibrational and rotational 
levels of impurities, eetees this treatment illus- 
trates that there is a lower limit on the resolving 
power of a tunneling spectrometer using normal metals. 
It is interesting to study the behavior of the 
total area or intensity of the second derivative of 


Gurrent vs joltage curves. “Fhus 
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The second integral in Eqn. (4.18) can be done exactly 


and we find 
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the second term in the integrand above is very small 
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compared to unity and we obtain 
Le.9) co 


2 
d 
| ao (SE) JdV =e | dw G(w) = constant . (4.20) 
0 dv 
0 
Therefore for spectra located far above thermal ener- 
gies the intensity or total area of the second deriva- 
tive of current should be constant as a function of 
temperature. This was also verified experimentally 


by Lambe and Jacklevic in the same paper. 


4.2 Inversion for G(w) in Terms of the Second 


Derivative of Inelastic Tunneling Current 


Let us consider Eqn. (4.14) for the moment as an 
integral equation for G(w) in terms of the second deri- 


vative of inelastic current which we shall denote by 
il 1 
TE (Vv) = ae (Giy” es (4222) 


For a fixed temperature T we shall assume that I"(V) is 
absolutely integrable, allowing us to define its Fourier 


transform by 
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ie ee Ai ' 
a Sketl ie een os | G(w) sin(sw)dw , (4.23) 
B V27 6 
where 
Kis) == | cos (Sx) Ki{x) “ax=". (294) 
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Making use of the identity 


Sin (ax) 608 (P%) ax a 4 aah )n (4.25) 
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we obtain 


= | eRAEy sin(ws)ds . (4.26) 
TT 


0 


K(s) may be evaluated in close form to yield 
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Therefore the barrier function G(w) may be determined 

from I"(V) by means of the formula 

es ines 
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The zero temperature limit is easily obtained making 


use of Eqn. (4.25) and we note that 
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in accordance with Eqn. (4.12). Eqn. (4.28) represents 
a valid inversion for the spectral quantity G(w) in 
principle. In practice the experimental noise and 
background inherent in measurements of I"(V) would 
severely limit its applicability. However, given the 
measured quantity I"(V,T') at some temperature T', it 
allows us to deduce the corresponding I"(V,T) at some 
higher temperature T >T'. If we substitute Eqn. (4.28) 
into Eqn. (4.19) a long calculation yields the follow- 
ing relationship 
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where the modified kernel function K(x;y) is defined by 


K (xpy)= 2 (e**(x-xcosy-ysiny) + e*(x+tysiny-xcosy) 

sy)=— 

y? es ~ 2e*cosy + 1) (e*- 1) | 
(4.30) 

In the zero temperature limit this reduces to the fami- 


Liar orm; 
RUS O). = K(x), (previously defined in Eqn. (4.15)). 


It may be noted that for certain values of the yatio 
T'/T the kernel K(x;y) is quite simple. For example, 


if T'/T = 1/2 then 
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ellen ee se 2 (23%) 
as (ey “e-rels) te Sani he ox 


K(x,1) corresponds to a distribution of height Be and 
linewidth 4.4k,T. Eqn. (4.29) may be of use to workers 
in low temperature tunneling spectroscopy who wish to 
isolate the natural linewidth broadening in sharp spec- 
tra due to the thermal smearing of the Fermi levels in 
the electrodes from other sources of broadening inherent 


in the barrier. 
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CHAPTER V 


TUNNELING IN SUPERCONDUCTORS 


5.1 The Semiconductor Picture 


The subject of electron tunneling can be dis- 
cussed in varying degrees of sophistication. The 
Simplest method is that used by Giaever (1960) to 
explain his original experiments in superconductive 
tunneling. In this method the superconductor is 
treated in terms of single particle excitations only 
(quasi-particles) since they are the ones involved in 
the tunneling process. The quasi-particle density of 
states is drawn with an energy gap of 2A centered at 
the Ferm: energy as "in Figs (57 la) 2° "At-T ="0°R -all 
states below the gap are filled and those above the 
gap are empty. At finite temperatures particle-hole 
pairs are created, the particles existing above the gap 
while the holes lie below the gap. Tunneling into the 
gap region is not allowed and this leads to structure 
in the I-V curve which gives a measure of the energy 
gap as well as of the quasi-particle density of states. 
In this approach to superconductive tunneling the 
superconductor is treated much like a semiconductor 


and consequently it is called the semiconductor model 
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of a superconductor. Giaever was able to obtain a 
good fit to experimental data by using this simple 
model. 

Following Giaever, we make the following assump- 


tions in deriving the tunneling current. 


1) The quasi-particles of both superconductors are 
independent particles satisfying Fermi-Dirac statistics 


and have an occupation probability given by 


free eT (Sack) 


where € 1S measured from the Fermi energy. 


2) The difference in Fermi energies between the two 


Superconductors is given by the applied voltage V. 


3) There exists a matrix element Mt connecting the 
states k and k' on opposite sides of the barrier. This 
matrix element is considered constant when the energy 


Use vanled initthe regiony offjinterest. 


4) The density of quasi-particle states is considered 


to be constant when the metal is in its normal state 
p; (e) = Py (5k 2) 


and is given by 
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py, (€) = 0, n(e) CB is) 


when the metal is superconducting. n(e) is the reduced 


density of states and according to BCS is given by 


| 


n(e) : Ven tame (5.4) 


= 0 Retpece i (Sos) 


where A is the energy gap considered to be independent 
of energy but not temperature. 

With these assumptions we can write down the ex- 
pression for the tunneling current in terms of tran- 
Sition probabilities in exactly the same way as we 
obtained Eqn. (3.39). 


(Ee yan nae (5.6) 


eae way 
In section 2.5 the matrix element was evaluated for the 
case of normal electrodes and it was found to contain 
density of states factors that cancelled the densities 

of states in the summations over right and left hand 
states. in superconductors it can be argued (Bardeen 
(1961) and Harrison (1961).) that the density of states 
fACtOrs an [Mal are the same as in the normal metal. 
In Sreconate ae tunneling it is necessary to make this 
assumption to obtain agreement with experiment. Thus 


we simply take 
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= dicnacty els! | (50) 


and Eqn. (5.6) becomes 


cl! 


tad | dE, [£° (Ey py) -£° (BE, -y+v) ] ||? 
S20 (ae si 
(ie 8) 


If we now make the replacements as in Eqn. (5.3) we 


find 
I=, fae [f° (ec) - £2 (c+V) Jn (e)n' (e+v) (5.9) 
where 

oO a 4te ) |m| ' = 8) 

=i 

In the limit that both metals are normal n(e) = n'‘' (e€+V)=1 
and 

A Oe a 


Therefore oe ise puSst tie conductivity of (the junction 
in the normal state. Using the definitions of the 


densities of states given by Eqn. (5.4) we obtain 
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where (2,r) denotes (left, right). We shall not pursue 
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the quantitative aspects of Eqn. (5.11) as these 

are given in a number of excellent reviews (e.g. 

see Douglass and Falicov (1964) or Solymar (1972)). 
Rather, we simply present the results in diagramatical 
fashion and discuss the results qualitatively. 

Fig. (5.la) shows a semiconductor type of ener- 
gy diagram for the case of tunneling between a normal 
metal and a superconductor. The density of states 
functions, according to BCS; are also included. At 
T = 0°K there are no excited quasi-particles in the 
superconductor and no excited electrons in the normal 
metal. Under these conditions it is easy to see that 
there is no current until a voltage of V = +A(0) is 
reached as in Fig. (5.l1b). At V = +A(0) the current 
rises rapidly due Caveat sharp Fermi energy in the 
normal metal and also the very high density of states 
in the superconductor at the gap edge. As V is in- 
creased further the current continues to rise and 
approaches the curve for normal-normal (N-N) tunneling, 
coinciding with it stor Vs> A(0)e This: 18. due to the 
fact that the superconducting density of states returns 
to that of the normal metal for energies e«>>A(0) (Eqn. 
ore 3) ae 

For finite temperatures there are excited quasi- 


particles above the gap as well as quasi-holes below 
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Fig. (5.1). (a) Energy level diagram for tunneling 


between a superconductor and a normal metal. (b) I-V 


curves expected on the basis of Simple theory. 
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the gap in the superconductor. In the normal metal 

the Fermi level becomes smeared. Thus when V is 

increased from zero, current immediately begins to 

flow and continues to rise gradually until V = +A(T) 

is reached. At that» point the density of states in 

the superconductor becomes large and the current 

increases more rapidly and approaches De FOr 22 AKt) 
The tunneling between two superconductors 

(A. # Ay) ns depueted in Fig. S(a2)scACtb 0°K there 

are no excited states and tunneling does not occur 

until a voltage of V = +(A, (0) +A, (0) ) is reached. 

At this point the current rises discontinuously because 

the density of states on both sides of the barrier is 

infinite at the gap edge. For V >> (A, (0) +A, (0)) the 

current approaches the n-n curve for tunneling as before. 
At finite temperatures and for V < (AL (T)-A, (T)) 

the current increases gradually because the number of 

quasi-particles available for tunneling is increasing. 

However in the range (A (T)—&, UE) < V < (AL (T) +A, (T)), the 

number of quasi-particles available for tunneling re- 

mains constant while, at the same time, the density of 

states in the right decreases. This leads to a decrea- 

sing current and thus negative resistance as shown in 

Pige t(o92b ie According: toikqn: (Susl bb) gat V=A.(T)-A, (T) 
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between two different superconductors. 


(b) I-V curves 


expected on the basis of simple theory. 
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V= A, (T) +A, (T) we again have an infinite density of 
states on both sides of the barrier and there is a 
discontinuous rise in current even for finite tem- 
peratures. 

The special case of tunneling between identical 
Superconductors is shown in Fig. (543) sm.Ac TT =0°% 
there 1s:no current Until Vi = 24 (0) § is reached, 
whereupon because of the infinite density of states 
on both sides, there is a discontinuous rise in current. 
At finite temperatures, as the bias is slowly increased 
the current slowly increases until the point where there 
can be no back current flow from the lower biased elec- 
trode. At this point the current remains nearly cons- 
tant. At low reduced temperatures (T/T, <<1) a-siaght 
decrease in current is expected because the quasi- 
particles from the left face a slowly decreasing density 
of states on the right. 

In real tunnel junctions the discontinuous current 
sump iat V = 2A4(T) has never been jobserved. Ihere are 
several possible reasons normally given why perfect 
agreement should not be expected. First, the prediction 
is based on BCS theory and it is unlikely that an ideal 
BCS superconductor really exists. Second, most of the 
reported experiments have been done on thin films where 


tunneling occurs into many randomly oriented crystallites, 
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giving a gap that is probably averaged over many direc- 
tions. For an anisotropic superconductor this would 
result in a smeared current rise. Another likely cause 
for this is the smearing of the gap due to a non-uniform 
strain across the tunnel junction. In the next section 
we show that this smearing may be partially explained 


by life time effects. 


5.2 Non-Eguilibrium Tunneling in Superconductors 


In Chapter III, we stowed! thet the zero bias 
anomalies in conductance seen experimentally were the 
result of blocking effects due to the finite relaxation 
times of the electrons. It is then natural to ask what 
effect finite relaxation times have on the tunneling 
Gf quasi-particies in superconductors. To calculate 
these effects, we use the semiconductor picture of a 
superconductor and proceed according to the assumptions 
in Section 5.1. However, we now write the expression 


for the tunneling current in the form 


Coie) 


where it is understood that the distribution functions 
in the above equations are steady state rather than the 


equilibrium ones. We may now proceed to define effec- 


pun a sso sient — ee 
cic ses ae ame res 
aoitses s¥en oid at sfighvsnet ‘aams acs ur en 
beaisiqzs ylisitieg sd ysm paizaome aint send a acl 
vesvarte mis REL Ye 


seid ores sit dens bowote ov . ILS wesqeni> at 
ait exsw ylistmemitsqxs ags2 scqastoubaoo ai esilamons 
Hodtsxsies stimtt ofd.os sub atostte piideeld to tives: _ 
Seilw des o¢ isavgen med? si 31 .anottosle one Io eamsz 


piifeangs sid ne over BSUS ‘cokjsxsiox o¢ina? Jtootite 


| ste ivolso oT puatliinid wigite ai aeloisisgq-iaseup to i, 7 
6 20 610d9hg zodDuUbAODImas Ss sey SW ,ag09326 saett : 7 
i ha beets oie Sit) oF ‘Pntbraens baaperg bas tosoubRoszequ a - 
noiaestaxs | sit stisvw won sw \tevewoH «1. 2 aoft0se a - 
orrek of at daetxup pakignay3 of? io? - ae 


(Stee) 


Ly. 


tive relaxation times for the quasi-particles in both 
Superconductors and derive the deviations of the 

steady state distribution functions from equilibrium 

as in Section 3.1. However, in that section the de- 
viations from equilibrium were derived for a system 

of independent particles (electrons) satisfying Fermi 
statistics and having in general an unspecified density 
of states. This is in accordance with our assumptions 
about the behavior of quasi-particles and we can simply 


take over the results of Eqn. (3.40) 


fr-f, 
phy a fy -fh en (Onno) 
aa sear = |M| (Oy, T * Px tt) 
where Ty and p, are the effective relaxation times and 


densities of state respectively for the quasi-particles 
in a superconductor. The summations in Eqn. (5.12) are 
then done in the standard fashion 
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We assume the superconducting quasi-particle relaxation 
times Ty depend only upon their excitation energies 


above the Fermi surfaces in the respective metals. 


Furthermore the densities of state are as given in 
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Eqn. (5.2). The integral in Eqn. (5.14) may then be 
replaced by one over the range of excitation energies 


€ and we obtain 


eee | de [£°(e) - £°(e+V) In(e)_n' (e+v) 
ad i +o [ml ? (oy mle) Ul (crv) om lerv it (6) 


== CO 


Cbs LS) 


where the primes refer to different sides of the junc- 
tion and Im|? and p, are the averaged matrix element 

and normal density of states at the Fermi surface. We 
note that if the relaxation times are zero then Eqn. 
(5.15) reduces to the conventional superconducting 
tunhnel-Ccurrent given by Eqni( 509). *Lakewise at both 
sides of the junction are normal, the reduced densities 
of states are unity and we obtain the normal ZBA current 
as in Egn. (3.42). We shall now proceed to discuss the 
superconducting relaxation times to be used in Egn. 
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543 Superconducting Relaxation Times 


Bardeen, Richayzen and Tewordt (1959) in the 
course of an investigation on the thermal conductivity 
of superconductors calculated the life times of quasi- 
particles due to elastic scattering processes. They 


found that the elastic scattering relaxation time of 
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a quasi-particle of excitation energy € above the Fermi 


surface in a superconductor is given by 


ar = sed elles oh Fy ea) 


where re is the relaxation time in the superconductor 
while Ty is the relaxation time in the normal metal. 
For € >> A the quasi-particles behave essentially like 
normal electrons in the sense that they have the same 
scattering times. However near the gap edge (ce = A) 
the life times are infinite. This is not as serious as 
it sounds, for in the same paper, they showed that the 
speed or group velocity of the excitations or quasi- 


particles ve is given by 


Las ae ee Cardo) 


where v,, is the Fermi velocity of the normal electrons. 


N 
Thus thesmean free paths of the quasi-particles due to 


elastic scattering is given by 
= = Ses 
QL ey, TV ( ) 


so that the mean free paths of the quasi-particles are 
the same as those of the normal electrons. Nevertheless, 


the fact that the life time is infinite is a measure of the 
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correlation the quasi-particle feels due to the pre- 
sence of the condensate although the elastic scatter- 
ing processes do not allow it to lose energy and thus 
recombine into the ground state. 

In the case of a quasi-particle in the presence 
of a lattice or phonon gas, Tewordt (1962) showed that 
the relaxation time of a quasi-particle of excitation 


energy € above the Fermi energy is given by the lengthy 


formula: 
eee mire a? 
meme || ese 
Semen 
le ee Gee a 
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yee pale eS 
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Physically, the first term corresponds to a quasi- 
particle state scattering into another quasi-particle 
state with emission of a phonon of energy w. The second 
term is the same process except that a phonon of energy 


w is absorped. The third term is the recombination 
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term and corresponds to the destruction of a pair of 
quasi-particle states of opposite momenta and spin 

with the creation of a phonon of energy w. It may 

be easily shown that in the limit of a normal metal 

(A = 0) Eqn. (5.19) reduces to the expression fon the 
life time obtained for a. normal.electron in Eqn. (3.26). 
The quasi-particle-phonon relaxation time given by 

Eqn. (5.19) is not infinite at the gap edge as in the 
case of elastic scattering because of the presence of 


the recombination time. 


5.4. . Phe Nonequilibrium Tunne lang wourrenk 


In order to calculate the nonequilibrium tunneling 
current we take our effective quasi-particle relaxation 


time t(e) to be given by 


where tT; is the normal metal elastic scattering time 
and the factor n(e) comes from the definition of its 
superconducting counterpart in Eqn. (5.16). T ose) Ls 
the gquasi-particle-phonon relaxation time. If we 
substitute Eqn. (5.20) into the expression for the 


superconducting tunnel current given by Eq. (5.15) we 


obtain 


isteom L[smion & to ae sit aL sana bal tiaras . 
att “or noLsaergee edd of e6ouhox (21.2) .apa (0 =a) : 
_(3¢.£) Lapa wk Aeaga6hs Lamson & toY Benkeddo emis out 
vd nS6vip omia moisaxsioxz avidodq-siol+asgq-iessp ot 


oid ui es sphe gee ods ts Stinktai fom et fete)” .np@ | 
to sonsse1q ond to Sevsced piitasdtave bitenle %o 9459 : 


4ints woitentdmeoes. sit ls” 


¥ om f 
_ 


" AeciD pail ennut dm ered topenott edt @u@ 6 
pritoanavgs me budid topecon afiy atetuoibs oF ~ebi0 at | Aa 7 
 noftsxslex afsttivsq~tesep svisoette 2vo Si6s aw smexa 
; . Vy 
ya mevip od of (3) said rn 
ies, tele | i 1 
act OA. Oy Br. ~ = 
qe : 
wale 
a 


amit pnistotispa olsesis ‘Istem Lamon. std at i orbdw . 
edi Yo noktigt29b sit} mos? sbmoo (a)n Zojos) ait Bas 9 
ai (3) gt (85/2) apa no: qxed au 3 “phitoubnoosequa 74 
ow 2t «pints dt sens Lox nonorg~stoidesg- ianinp fasup oid = 
Pascale AOR as st aps sutiseaia ae: 


wa. 


ee 


hia? de [f° (e) = £°(e+v)).n(e)n' (e+Vv) 
oO t ' 
me ‘Earn Zan . 
ee et He 4, AE eB rs 
‘deal e+V) eee. 
pl cate eke. at ie ae 
Cops? eot 


Uae) 


where the primes refer to the opposite side of the 


junction and tT, is defined by 


B 


1 27 Mamie: 
ee ah ele 
aS in Eqn. (3.47). We may note that the tunneling 


current given above is an odd function of bias since 

T ap ‘€) is an even function of its argument. Eqn. (5.21) 
may be written in a simpler way by assuming that we have 
identical superconductors on both sides of the barrier. 


In this case 


£P(e) - £° (e+v) 
piney moonlit apes 8 oy en 
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(52.2) 


The above expression for the superconducting 
tunneling current exhibits a new feature which is not 


present in the conventional current given by Eqn. (5.9). 
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That is, the tunneling current does not have a dis- 
continuous jump at V = 2A(T). Rather, if the ratio 
tTi/Tp rs small, tt starts increasing at Vi =—200)) wise 
a large but finite slope. For V >> 2A(T), we know trom 
the behavior of the ZBA in normal tunneling that the 
factor in the denominator above is very small and hence 
the tunneling current reverts back to the conventional 
BCS current. The reason we do not get a discontinuous 
jump in the tunneling current at V = 2A(T) is due to the 
fact that the life time factor in the denominator of 
Eqn. (5.22) can never be zero. Therefore the infinite 
density of states at the gap edge which causes the jump 
mec urerente LSanO tet Cilkt. 

We shall now consider the case in which the elec- 
tron-phonon interaction is weak (i.e. the recombination 
time is large compared to the elastic scattering time). 


Thus the term given by 
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is of the order of unity and the tunneling current takes 


the form 
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The T = 0°K tunneling current Wis*plotted in Fig. (5.4) 
for various values of the ratio T,/T. and compared to 
the conventional BCS tunneling current. We may remark 
that: the £1 for the value of Ti/Tp x .001 corresponds 
Closely to the experiments in superconducting Al-Al 

junctions by Giaever, Hart and Megerle (1962). Assum- 
ing that the nominal resistance of their junctions is 
of the order of 502 so that the tunneling time TR Can 

ul 


be taken as tT, * 10 ° sec., we obtain-the value 


t, x 10°79 seconds . (5.24) 
This life time agrees with our ZBA results in Chapter 
III. It is interesting to note that in the same paper 
Giaever et al attempted to account for the finite slope 
of the current rise by assuming that the energy gap was 
smeared. They find that if they attribute the smear- 


ing to, a Ife time, they obtain a) lower Jimit or tome 


seconds. 

The smearing of the current rise as given by 
Bon. (5.22) or Eqn. (5.23) depends upon the magnitude 
orether ratio T;/Tp- Since Tp is a measure of the re- 
sistance of the junction we disregard it in the follow- 
ing discussion. The elastic scattering time v should 


depend strongly on the properties of the metallic films. 


In practice we expect the amount of elastic scattering 
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Pig. (5,4). Tunneling current due to life time effects. 
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to increase significantly as we either decrease the 
film thickness of the electrodes or else increase 

the impurity concentration of the films. Thus Ts/T. 
should decrease accordingly. The width of the current 
rise should then decrease as the film thickness is 
decreased and the impurity concentration is increased. 
This is in agreement with the observations of Rowell 
(1969) who found that the width of the current rise in 
Pb junctions decreased as the film thickness was reduced 
from 2000 to 1000 A, and the current rise became very 
sharp when In was alloyed into the Pb. 

In the case of finite temperatures, the calcula- 
tion of the tunneling current is numerically much more 
complicated when we take into account the relaxation 
time due to the lattice and the temperature dependent 


gap. ‘This work 1s still in progress at this time. 
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CHAPTER VI 


GEOMETRICAL EFFECTS IN TUNNEL JUNCTIONS 


AND CONCLUSION 


6.1 Structure of the Tunnel Junction and Splitting of 
the ZBA 


In previous chapters we treated the tunnel junc- 
tion under the assumption that it was a homogeneous 
sandwich of the two electrodes and the insulator. In 
reality, tunneling occurs only ateecertain portions of 
the cross strip area of the junction. This was conclu- 
Sively shown by Adler and Kreuzer (1972). They alter- 
nately passed a large current in the top film (Tt ? and 
the bottom film of therjunction (I) while simultan- 
eously measuring the tunnel current (i,) and its 
derivatives (Fig. (6.2a)).. “Fig. (6.2)s shows a typical 
set of conductance curves that they obtained in an 
Al-Al junction at 4.2°K. When the conductance was 
measured in a conventional fashion without any addi- 
tional currents flowing in either the top or bottom 
fiim-the solvda-curve-in Fig. (6-1) "was Obtained. 

The ZBA dip in conductance is seen very clearly. When 
a large current (T.) was passed along the top filn, 
the ZBA split into two minima as shown in the dashed 


curve. However, upon passing a large current (I) along 
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Fig. (6.1). Conductance as a EUNnCcbLon of bias ior an 
Al-Al junction at 4.2°K: (——, no current through 
Seitherrelectrode) ;ate =, ;5SeAvsalongethe top, film) ; 


(<=, =. =, .5 A along the bottom film). 
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the bottom film, it was found that the ZBA character- 
istic smeared out as shown by the dotted curve. 

This experimental result may be easily explained 
if one assumes that tunneling occurs only along the 
edges of the bottom film. Thus when the current I 
is passed along the top film, there is a potential 
drop along the film and the two effective areas for 
tunneling differ in potential. This causes the ZBA 
characteristics to now have two minima whose centers 
are separated by the potential drop along the top film. 
On the other hand when a current I, is passed along the 
bottom film, there is a continuous potential drop along 
the edges and the ZBA becomes smeared out. Adler and 


Kreuzer summarized their results as follows: 


uk) Splitting? of jhe~ZBAsas ondysopseuveduwrehaa 
current : (this has been tested with currents in 
excess of .5 A). Similar results have been obtained 
using Ag,alny, Sm, andarbeclectrodes: 

2) In all cases the current 1, splits the ZBA mini- 
mum into exactly 2 minima. 

a) The shift of the minima to either side is always 
a linear function of Iv. This ohmic behavior was found 
to occur over the whole range of metal film thickness 
used (200-4000 A). Such ohmic behavior was also tested 


in several triple junctions having a common bottom 
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film and oxide while the top films varied in thickness 
DY alratio or 122.4. 

The above results clearly suggest that the oxide 
barrier has two thin regions along opposite sides of 
the bottom film. The fact that no splitting arises 
from the current I, in the bottom film then indicates 
a rather uniform barrier parallel to the bottom film. 
The junction may then be idealized by the geometry 
shown in Fig. (6.2b). The presence of a current I, 
along the top film may then be represented by two thin 
regions at y = tb being biased at two different poten- 
tials tV, as in Fig. (6.2c). The magnitude of the 
potential difference ve is then determined by the 
current I. and the film resistance. When a current I. 
is passed along the bottom film, there would be a con- 
tinuous potential drop along the constant thickness 
tunneling region as indicated in Fig. (6.2d). Adler 
and Kreuzer assumed that the tunneling conductance per 


unit area of the thin regions was given by 


o(V) = oe + 60(V) 


where 60(V) is the negative contribution to the con- 
ductance due to the blocking effect and is given by 
Egn. (3.52). For simplicity, they considered 66 (V) 


to be given by the T = 0°K limit as in Eqn. (3.56) 
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Schematic of the junction geometry as described in the text. 
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and showed that both the smearing and splittings of 
the ZBA characteristic could be accounted for by 
uSing the idealized geometry in FG (O62). 

It follows that tunneling occurs predominantly 
near the edges of the junction and we are justified 
in assuming that the effective area for tunneling is 
a small fraction of the cross strip area as indicated 
in evaluating Eqn. (2.22). The method of introducing 
potential drops across the metal films by external 
currents may be used to probe the structure of the 
oxide. In principle it should be possible to obtain 
the conductance per unit area as a function of posi- 
tion in the junction by carrying out an exhaustive 
set of measurements with varying I, and Iy: By using 
superconductor in the bottom film the resolution of 
such measurements with I, may be improved by using the 


gap edge as a potential probe rather than the ZBA in 


the normal metal. 


6.2 Conclusion and Discussion 


In normal metals, we have shown that the observed 
structure near zero bias in the conductance o and its 
derivative do/dV can be accounted for by including the 


effects of the finite electron relaxation rates in the 
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junction electrodes. The non-zero life times then 
lead to an observable blocking of otherwise available 
electron tunneling states at low bias. 

The dependence of the conductance minima on the 
properties of the junction may conveniently be sum- 
marized with reference to the T = 0°K diagrams in 
Figs. (3.3) and (3.4) for the models n=l and n=2 res- 
pectively. We note that the amplitude of the conduc- 
tance minimum is proportional to the ratio t;/Tp- 

Hence by increasing the amount of elastic scattering 

in the electrodes we would decrease this ratio. 
Similarly, making the oxide thicker would increase the 
tunneling time Tp and again make this ratio smaller. 

A possible experiment to test the dependence of the 
size of the conductance dip would be to construct a 
multiple junction with the same base metal film and 
oxide. The films deposited on top could then be alloyed 
with a suitable metal in increasing concentration. If 
our conjecture is correct, then the size of the conduc- 
tance dip should decrease with increasing concentration 
Of the alloy material in the top f1im. » To (test the 
dependence on oxide thickness would be more difficult. 
However, by carefully examining the structure of the 
oxide as described in Section 6.1 one should be able 


to, in principle, also test this aspect of the theory. 


a 4 ie va a! 
nowy aqmit stil orss—on otit’ 


aldalisvs Ssatwiedto 36 ibtsot nhigld ofl 
>> - 
anid ead te getese ened: 


eis ao smindna sonstoubaos | eit to soasbneqsh * 1 
~mia od yitneinevnbo yem notgono[ ot io eels oa 
ni dps ib HO = T edt of Sclinuwinte ut iw ovis ee 


-@9% Son bas f=n aisbom ei3 x09 eel saiee (E. é) me i 


7 
~oubaor Sd3 to sbugifqms sit tsdy ston a /qlevisoeq a 


: 
at Net ofjet sat oF iene init ee li ai moninim sor6t if. - 
> ae 
r] 


puristsss2 oitjesis Yo tavoms ors piiasstont vd soreH 
.olisx sti} oecoiseb Bivow sw eeheiszeis ont ak 

end Sessz5ni £ Luow soileids sbixe old pattem yie tim]: 
»weilsme oitst eins sdsm tisps bis 5) smist onktennws 


edt Zo sunshnaqeb ont tesy 09 Sasatsequs aftieaeq A - 
« 


# Jourdane oF ed Biow Web Lah pibsranibends std io oala 


Sms mist ledom Sasd omée oo rigiw Aolioas, eiqitiom 
fie vests 43 nods bivon qo? no Setiecges aati? aT - 8bLx0. 
pat -noitpsineonos agi ai Ladom etdas tue 6 ow -_ 
“—ogbneo sit to sgie od cata f292200 ai exwdoetaoo 300: a 
noitsatnoeonop prizsoxomt ris iw vemee bivode ath ones a 


ats. daot of “ym 2 abe eae ai isizosae yolts ott to 
_ 
-tinoittib stom Sd sivow, aeonitoins ehixe no eonebas a9b 


& 7 


git 2o ora ede, erld parte “s yiluieres yd xa ws Bf 
-oid& ed biuoda sno Le a wai 9302 at bodizp2eb ae 4 ad ‘ 


edt to tosqes. ei Sil ois at ates f t Ces 
iL? a - oa? ¥ a in 


: yxoant 
. 


134 


The width of the conductance dip as indicated in 
Figs. (3.3) and (3.4) is proportional to the quantity 
Ce ie where n=l or 2 according to which model 
we choose for the electron-phonon coupling. Thus as 
we increase the coupling strength of the electrodes 
or else the elastic scattering time i the width of 
the conductance dip should decrease accordingly. A 
much more difficult experiment to perform would be to 
construct a tunnel junction with a crystal on one side. 
In this case Ts should be very large compared to its 
value in the evaporated film and one should be able to 
observe a fairly deep and sharp conductance dip. Also 
by tunneling into single crystals one should be able to 
probe, in principle, the anisotropy of the electron- 
phonon relaxation times. These considerations can like- 
wise be carried through for the derivative of conductance. 
The agreement of theory with experiment seems to 
indicate that the model aml nDNG e « w iS a more rea- 
sonable choice to make. However a better test to dif- 
ferentiate between the two models should be possible 
at low temperatures. Fig. (3.4) indicates that for 
a) ee x a there is a kink in the derivative of 
conductance at zero bias and zero temperature. This 
feature rapidly disappears at higher temperatures but 


calculations indicate that it should become noticeable 


below about .5°K. This may serve as a sensitive test 
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for the two models although modulation problems may 
become severe in this low bias regime. We have no 
physical explanation for this kink except to say that 
iti ais a resultsof thesfactithathdg/dVegoes like v* at 
very low bias and at zero temperature. Although we 
have not done so, a more complicated version of 

ciated alucl} may probably exist to fit tie theory ina 
more suitable fashion. However, there seems to be 
little precedent in present unrelated theoretical or 
experimental work tO_ JUStALY this. 

In the work we have done on superconductive 
tunneling, the nonequilibrium tunneling of quasi- 
particles seems to indicate that the tunneling current 
between two identical superconductors does not have a 
discontinuous rise. This is in agreement with experi- 
ment. However, we have used the semiconductor picture 
of a superconductor both in the derivation of the 
transport equations for nonequilibrium tunneling and 
the tunneling current itself. A more rigorous treat- 
ment using the microscopic many-body theory of super- 
conductive tunneling in the presence of relaxation 
processes would be in order. Nevertheless, it is en- 
couraging that the same choice of life time parameters 
used to explain the normal ZBA also serve to give quali- 


tative agreement with experiment in the case of super- 


conductors. 
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